Hamiltonian theory of gaps, masses and polarization in quantum Hall states: full 

disclosure. 
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In two short papers I had described an extension, to all length scales, of the hamiltonian theory 
of composite fermions (CF) that Murthy and I had developed for the infrared, and applied it to 
compute finite temperature quantities for quantum Hall fractions. I furnish details of the extended 
theory and apply it to Jain fractions v = p/{2ps + 1). The explicit operator description in terms of 
the CF allows one to answer quantitative and qualitative issues, some of which cannot even be posed 
otherwise. I compute activation gaps for several potentials, exhibit their particle hole symmetry, the 
profiles of charge density in states with a quasiparticles or hole, (all in closed form) and compare to 
results from trial wavefunctions and exact diagonalization. The Hartree-Fock approximation is used 
since much of the nonperturbative physics is built in at tree level. I compare the gaps to experiment 
and comment on the rough equality of normalized masses near half and quarter filling. I compute 
the critical fields at which the Hall system will jump from one quantized value of polarization to 
another, and the polarization and relaxation rates for half filling as a function of temperature and 
propose a Korringa like law. After providing some plausibility arguments, I explore the possibility 
of describing several magnetic phenomena in dirty systems with an effective potential, by extracting 
a free parameter describing the potential from one data point and then using it to predict all the 
others from that sample. This works to the accuracy typical of this theory (10 -20 percent). I 
explain why the CF behaves like free particle in some magnetic experiments when it is not, what 
exactly the CF is made of, what one means by its dipole moment, and how the comparison of theory 
to experiment must be modified to fit the peculiarities of the quantized Hall problem. 
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I. INTRODUCTION 

There is a consensus among theorists and experimen- 
taUsts that the concept of the Composite Fermion (CF)lil 
is very useful in uruierstanding the fractional quantum 
Hall effect (FQHE)B. This concept allows us, among 
other things, to decide which fractions are robust (the 
Jain series), to generate very accurate trial wavefunc- 
tions and gaps, and to determine the allowed values of 
polarization when the spin is not fully polarized. 

The aim of this paper is to describe, in detail, a hamil- 
tonian formulation which provides a comprehensive way 
to describe FQHE states qualitatively and quantitatively, 
at zero and nonzero temperatures. Recall that in the 
theory of superconductivity or Fermi liquids one always 
seeks a transformation relating the original electronic 
variables to those of the ultimate quasiparticles, these 
being the Cooper pairs and Landau's quasiparticles re- 
spectively. In the hamiltonian formalism used here one 
passes form a description in terms of electrons to one 
in terms of CF's through a sequence of transformations. 
One ends with an operator description of the CF. Even 
though such a change of basis is approximate, it provides 
us with valuable insights, and occasionally, quantitative 
information that is unavailable in the wavefunction ap- 
proach, such as ways of coupling to impurities and cal- 
culating unequal-time correlations. 

Given that the FQHE has no small parameters, how 
was this passage possible, even approximately? The an- 



swer has two parts. Originally Murthy and M developed 
the transformation using a combination of the random 
Phase Approximation (RPA) and the infrared limit, and 
obtained the electron density operator in the CF basis 
at small ql, q being the momentum and I the magnetic 
length. The RPA kept operators at different q^s from 
mixing during the transformation. Given the density op- 
erator, the hamiltonian, which is just the interaction en- 
ergy when electrons are restricted to the lowest Landau 
level (LLL), could be written down. For potentials that 
were soft atJarge ql, we could compute objects like acti- 
vation gapso. 

I then extended the infrared results to all ql by ap- 
pealing to certain algebraic properties known to be true 
for the LLL projected problem in any basis. The exten- 
sion consisted of taking the small ql series for charge and 
the constraints, and exponentiating them. The resulting 
operators obeyed the desired algebras. Given that the 
extensions were not unique, but just minimal and consis- 
tent, it was very satisfying that they embodied, despite 
their questionable pedigree, much of what was known 
about the internal structure of the CF, illustrating once 
again that the tight constraints of the FQHE problem 
in the rlttL can actually work in our favor. Two short 
papersB-Q described this extension and its application to 
finite temperatures, T > 0. Here I provide the promised 
details and additional insights gained in the meantime. 

The topics covered here fall into two classes. 

The first has to do with matters of principle. For ex- 



ample, can CF be free or nearly so? Without going into 
any details, we can say no, since it takes two very different 
masses nia and rn^, to describe their polarization and ac- 
tivation phenomena, something impossible in a free the- 
ory. Yet certain polarization phenomena at T = seem 
to be very accurately fit by free particles of mass nip. The 
hamiltonian formalism not only allows one to compute 
these distinct masses, it also resolves the paradox posed 
above. Next, CF are supposed to derive their kinetic 
energy from electron-electron interactions. The present 
formalism provides an explicit expression for not only this 
kinetic energy but also the CF interactions, i.e., the full 
CF hamiltonian. Both kinetic and interaction terms have 
unusual functional forms, which are determined uniquely 
by the theory. Those attempting to fit data to traditional 
forms of energy should bear this in mind. Having a con- 
crete hamiltonian also eliminates questions such as which 
of nia and rrip should be used at T > 0. A Hartree-Fock 
calculation on the full CF hamiltonian gives the answer. 

The second set of topics consists of application aimed 
at showing that it is possible to compute, by analytic 
means and often in closed form, numerous physical quan- 
tities pertaining to FQHE states. I consider here the 
computation of activation and polarization masses and 
gaps, charge profiles of quasiparticles and quasiholes, 
critical fields for magnetic transitions from one quantized 
value of polarization to another for the gapped states, all 
of which arc T = quantities. I also compute the po- 
larization P and relaxation rate 1/Ti at T > for the 
gapless fractions. 

Given the importance and utility of the CF idea, and 
the potential for misunderstandings, I have made every 
effort to make my arguments accessible to as wide an au- 
dience as possible, underscoring the various assumptions 
that go into the calculations, emphasizing not only the 
formalism but the physical picture that goes with it. 

In Section II, I describe how one arrives at the hamil- 
tonian in terms of CF, starting with electrons. This dis- 
cussion will be brief, given that details have already been 
published. It does include recent insights on the internal 
structure of the CF, such as what it is made of and what 
exactly its dipole moment means. In Section III, the fi- 
nal equations are analyzed to gain familiarity with their 
main properties. The theory is then recast in a form that 
makes it more amenable to the Hartree Fock approxima- 
tion. Readers not interested in looking under the hood, 
may begin with the equations listed towards the end of 
Section III, which form the basis for the subsequent cal- 
culations. 

In Section IV, I calculate the activation gaps Aa for 
several fully polarized fractions within the HF approxi- 
mation and cepipare to the results of Park, Meskini and 
Jain, ( PMJ)lI based on trial wavefunctions. With one 
exception, all calculatioiis will be carried out using the 
Zhang-Das Sarma (ZDSp potential 
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here as a free parameter. I use this potential to illustrate 
the method, which is instantly adaptable to any other. 

This present calculation differs from earlier worka 
based on the infrared theory in that it yields finite re- 
sults even for the coulomb case, A = 0. The numbers 
agree to within 10 — 20% (and occasionally better) for 
potentials that seem to describe real systems (A ~ 1 — 2). 
Similar results are found when I compare to the exact 
diagonalization results of Morf et au for a similar range 
of the parameter h that enters their potential: 
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where I is the magnetic length and ZA = A was originally 
introduced to describe sample thickness, but employed 



Why bother to reproduce numbers that are already 
known, to a lower accuracy? The point is that the present 
approach is fully analytic, makes the underlying physics 
very transparent, and furnishes an explicit operator de- 
scription of the final quasiparticles, which permits a pre- 
cise formulation of many question pertaining to them 
that would be otherwise nebulous. It is also worth men- 
tioning that the closed expressions for physical quantities 
allows them to be computed in a few seconds on a PC. 
The reduced precision is a price we must pay in return. 

I compare the theoretical activatioB-gaps to the ex- 
periments of Du et a&B and Pan et aO. In comparing 
the theory to experiments, one needs to decide how to 
handle disorder and LL mixing, which are suppressed in 
the PMJ and Morf at al computations. There does not 
exist at present an analytical theory for incorporating 
disorder. (There does exist numerical work deiaopstrat- 
ing the effect of disorder, see for example Ref£3). My 
approach has been to take experimental points and fit 
them to the theory with the ZDS potential and ask what 
A is needed. This is done solely to get a feeling for its size 
and also compare it to the values computed for the pure 
system with no LL mixing, using say, the Local Density 
Approximation (LDA)E3J3. It is not assumed that the 
ZDS potential actually describes the problem at hand. 
It simply taken as a reasonable variant of the coulomb 
potential with a free parameter that can parameterize 
sample thickness and illustrate the hamiltonian method. 
When results for polarization phenomena are compared 
to experiment, a more ambitious approach to A is under- 
taken. __ 

I explore the question raised by Pan et aO of how 
the normalized effective mass of CF near half filling com- 
pares with that near quarter filling. Is the rough equal- 
ity, observed experimentally, in accord with theory (in 
the absence of disorder)? In general the answers depend 
on how the fractions are reached- by varying the den- 
sity, the field or a combination of both. Typically these 
masses lie within a factor of two of each other and there 
appears to be no deep reason wlui,they should be exactly 
equal, a point also made in Ref.EJ 

I provide the profiles of charge density in some gapped 
states with one quasiparticle or quasihole and compare 
to the unpublished work of Park and Jain based on trial 
wavefunctions. I explore the v = 1/2 case, especially the 
dipole moment and what it means, in some detail. 



Since the CF hamiltonian naturally separates into a 
free part Hq and an interaction Hj, I explore the effect 
of turning off Hj and find it can change the answer by 
as much as a factor of two. 

I ask how well particle- hole symmetry works, i.e., to 
what extent gaps sXv= p/(2ps + 1) in the fully polarized 
case equal those at 1 — z/ and find it works very well. I 
point out that this was not a foregone conclusion since 
the formulae for the two cases are quite different. 

Section V is devotedpto spin physics at T = 0, an area 
investigated in the pasttj and more recently by Park and 
Jainij in the wavefunction approach. In the absence of 
an overwhelming Zeeman term, one has to consider CF of 
both spins. The polarization of the ground state will be 
decided by a competition between ferromagnetism and 
antiferromagnetism. When the energy difference (not 
counting the Zeeman energy) between two ground states 
of different polarizations equals the corresponding Zee- 
man energy difference, a transition will take place. The 
transition can be driven, for example, if the density and 
field are varied together at fixed filling fraction, or by tilt- 
ing the sample at fixed perpendicular field and density. 
The critical fields B'^ at which these transitions happen 
are calculated. The calculations reveal a feature noticed 
by Park and Jain using trial wavefunctions, namely that 
they may be fit very well to free fermions with a constant 
polarization gap, Ap. How do we reconcile this with the 
fact that activation gap A^ is substantially different from 
Ap , and that turning off Hi makes a sizeable difference to 
Aq? I will show how two-dimensionality and rotational 
invariance can conspire to mimic free-field behavior for 
these polarization phenomena. For example, in the gap- 
less case oiv— 1/2, I will show that while the HF ener- 
gies £(k±F) of fermions on top of the spin up/down Fermi 
seas are not even quadratic functions of the correspond- 
ing momenta k±p, (and have substantial fc^^ pieces), the 
energy cost of transferring a particle from the top of one 
sea to the top of the other (which is what determines the 
polarization) takes the free-field form {k'^p — k'^p)/2mp. 
A similar situation exists for the gapped fractions. These 
arguments should caution experimentalists and theorists 
against misinterpreting the free-field fit. 

My results for the critical fields. B^ are then compared 
to the data of Kukushkin et aO. In the case of mag- 
netic phenomena I take a different approach to A. First 
A at one transition is obtained by fitting to the observed 
-B"^. Scaling laws then determine it for the other transi- 
tions, whose B'^ can be predicted to within 20%. In other 
words, it seems to be true for magnetic transitions in the 
samples considered, that the disordered system can be 
described by an effective translationally invariant poten- 
tial. I describe a limit in which this result can be jus- 
tified. However, I pursue this approach for all magnetic 
phenomena, even though at present I cannot provide sim- 
ilar arguments for all of them. I do so because it works 
to within the accuracy typical of this theory, a feature 
that needs to be understood. 

Section VI considers magnetic phenomena of gapless 
states as a function of temperature T. The T > physics 
is the first instance the present methods outperform com- 



plementary approaches based on exact diagonalization 
(limited to small systems) or trial wavefunctions (limited 
the ground state and very low excitations.) It has no 
finite size effects since one works in the themodynamic 
limit all along. One need not agonize over whether nia or 
nip should be used in computing a T > quantity such as 
polarizations: given a concrete hamiltonian, a Hartree- 
Fock (HF) calculation gives the results. The HF works 
well because most of the right, nonperturbative physics 
is already built into the hamiltonian. 

A Hartree-Fock calculation gives the polarization P 
and relaxation rate 1/Ti as a function of T and the poten- 
tial. These numbers are then compared to experiment, 
again by assuming that the real system can be described 
by an effective potential, fitting A at any one data point 
from each sample and explaining the rest of the data from 
that sample at other fields, tilts and temperatures. Into 
these calculations go the noncanonical, nonconstant, den- 
sity of states peculiar to this hamiltonian. These pcesults 
are compared to experiments of Dementyev et a£3 who 
measured P and 1/Ti at zero and a 38.3 ° tilt for a range 
of temperatures. They had pointed out that attempts to 
fit all four graphs with standard a hamiltonian (with a 
mass TO and Stoner coupling J) led to four disjoint set 
of values. On the other hand if the hamiltonian theory, 
with its peculiar functional form for H is used, a good fit 
to all four graphs is possible with a single A — 1.75. 

^n comparing to the polarization data of Melinte et 
I find the predictions work for the untilted case but 
not at a tilt of 61 °. The reasons for this are discussed. 

I provide an approximate expression for 1/Ti as a func- 
tion of temperature. It has the Korringa form only in the 
critical case where the polarization saturates exactly at 
T = 0. 

Conclusions follow in Section VII and end with a dis- 
cussion of a procedure for comparison of theory to ex- 
periment that is tailor made for the Hall problem. Many 
details are relegated to the Appendix, which ends with 
a summary of symbols for which there does yet exist a 
uniform convention. 



II. THE HAMILTONIAN FORMALISM 

Let us begin by tracing the path from the hamiltonian 
in terms of electronic coordinates to that in terms of CF, 
focusing on the spin-polarized case for fractions of the 
form 



P 



2ps + l 



(3) 



The results extend easily to z^ = ^ ^_^ . The treatment 
of old pubfished steps, presented here for completeness, 
will be schematic. 

Consider electrons of band mass m and number den- 
sity n, described by the following first quantized hamil- 
tonian: 
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V X A = -eB 
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where h = c = I, z the unit vector along the z-axis, I 
is the magnetic length, B is the applied field, V is the 
inter-electron potential, and tj is the cyclotron coordinate, 
whose components are canonically conjugate: 

[Vx ,%] = *^^- (10) 

Thus the spectrum is given by Landau Levels (LL): 



E = ujo{n+l/2) 
ujQ = eB/m. 

In the lowest Landau level (LLL), 

(?7>-LiL = I- 



(11) 
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There is a huge degeneracy of each LL due to the fact 
that the guiding center coordinate 



R 
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r — r z X p 



whose components obey 

[Rx , Ry] 
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does not enter H. The conjugate pair [Rx,Ry) ranges 
over the entire sample, whose area is its phase space, 
and determines the degeneracy if one employs Bohr- 
Sommerfeld quantization with l^ playing the role of fi. 
At the Jain fractions, the inverse filling fraction 



eB 1 

= ^2s+- 

2'Kn p 

= flux quanta per electron 

= states in the LL per electron. 



(16) 
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If i^ < 1, there is enough room in the LLL to fit all 
the electrons in the noninteracting case. One expects 
that if the cyclotron energy loq is much larger than the 
interelectron potential, the ground state and low lying 
excitations will be formed out of states in the LLL. 

Since 



r = R + r; 

a natural projection to the LLL is 
7^ ^ r ^ R. 



(19) 
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After this projection the two commuting coordinates be- 
come canonically conjugate. 



Given the huge degeneracy of the LLL for v < \, the 
problem is the selection of a unique ground state. Laugh- 
lin blazed one trail, writing down inspired trial wavefunc- 
tions for v = l/(2s-f 1). The other route is to try to start 
with the electronic hamiltonian and try to reach, through 
a sequence of approximations, the final quasiparticles, 
which in this work, will be the composite fermions. 

For Laughlin fractions, where the wavefunction is ex- 
traordinarily compact and simple, one-has the option of 
working with Composite Bosons (CB)E!j, which have con- 
siderable appeal of their own. 



A. What is a Composite Fermion? 

What exactly is a CF composed of? I am grateful to 
G. Murthy for some very useful discussions of this issue. 

Laughlin showed (using arguments involving adiabatic 
introduction of a flux quantum) that ai v — 1/ {2s -I- 1) 
the elementary excitations have a charge 
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(2s -fl) 
Consider the following state 
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The prefactor Y\Az — zo), multiplying Laughlin's ground 
state wave function ( whose gaussian factor has been sup- 
pressed) is a vortex. Due to the zero at zqj there is a 
charge deficit near that point, whose value, in electronic 
units, can be shown to be — l/(2s -|- 1). 

In CF theory the quasiparticle is believed to be an elec- 
tron bound to 2s vortices. We shall see that while this 
is clearly so for the Laughlin series, the situation for the 
Jain series is more complex. 

If we look at the Laughlin wavefunction we see a 2s-|-l- 
fold zero at each electron: one is due to the Pauli prin- 
ciple and the other 2s represent the zeros due to the 
captured vortices. There is no question of which vortex 
belongs to which electron since the vortices are sitting 
on the electrons. The charge of the electron plus 2s— fold 
vortex, i.e., CF charge is given by 
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(23) 



Sometimes the vortex is incorrectly used interchange- 
ably with a flux tube, 



flux tube — 
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(24) 



which has the phase of the vortex but not the zero. In 
other words the CF is described as an electron bound 
to 2s flrp, tubes. It was emphasized very early on by 
HalperincZl that for Laughlin fractions, electrons like to 
bind to vortices due to the couloipb attraction. This was 
also at the heart of Read's worked, which extended the 



concept to i^ = 1/2 where the wavefunction is obtained 
by projection to the LLL (more on this shortly) which in 
turn causes the vortices to move off the electrons. In any 
event, electrons are not attracted to flux tubes, which are 
neither charged nor low energy excitations. 

Consider the Jain wavefunctions &i v = p/ {2ps +1): 
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The factor Xp(z,z) describes p-fiUed CF Landau levels 
and the Jastrow factor J {2s) = Wi^j {^i ~ ZjY^ describes 
vortices sitting at the locations of the particles and P, 
the LLL projector, replaces the z's as perO 
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(26) 



Let us first ignore V. Then there are indeed 2 s zeros 
per particle in the Jastrow factor (located on the elec- 
trons) and one (not necessarily analytic) zero per parti- 
cle in Xp{z,z). (By zeros of the wavefunction, I always 
mean as a function of one coordinate, all others being 
held fixed.) The 2s-fold vortex has a charge 
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{2ps + 1) 



(27) 



in electronic units, a result that can be deduced from just 
the Hall conductance and incompressibility of the state. 
If we add the vortex charge to that of the electron, we 
do indeed get 



1 



2ps + 1 ' 



(28) 



J, is the correct, quasiparticle 
ictions. This is also confirmed 



which according to S 
charge at all gapped frac 
by focusing experimentsEj. So we may say on the basis 
of this unprojected wavefunction that the CF is the union 
of an electron and a 2s-fold vortex. Since the vortices sit 
right on the electrons, there is no confusion on who they 
are bound to and all moments vanish except the total 
charge. 

A lot of this changes upon projection by T': ^ acting 
on the Jastrow factors, moves the zeros away from the 
particles and many of them vanish, leaving 2s -f l/p per 
electron, (determined by the number of flux quanta per 
electron). Thus, after projection, vortices cannot be as- 
sociated with the electrons in an unambiguous way. For 
example at 2/5, there are 2.5 zeros per electron. One 
sits on the electron due to the Pauli principle, leaving 
1.5 non-Pauli zeros per electron, which will neither lie on 
the other electrons, nor be numerous enough to form two 
vortices per electron. Presumably, in states involving 
projection, where the wavefunction has a very compli- 
cated form due to the action of V, there is some nontriv- 
ial sharing of vortices between electrons. In the limiting 
case oi v = 1/2 only one non-Pauli zero per electron will 
remain after projection. We shall return to i/ = 1/2 later. 

These remarks do not imply that the CF approach 
to writing down trial wavefunctions based on electrons 



binding to vortices is in jeopardy. Thinking in terms of 
vortices still gives the unprojected wavefunction. This is 
all one needs, since the act of projection, while compli- 
cated, has a definite algorithm that is routinely carried 
out, and yields wavefunctions with their incredible over- 
laps with exact results. My message is only that in the 
end, if one looks at the projected wavefunction, (which 
is going to very complicated) one is not likely to find any 
simple correlation between electrons to vortices. These 
remarks apply to any projected wavefunction '^ lll- 

It is quite remarkable that even though we cannot as- 
sign to each electron a 2s- fold vortex in ipjaim e* is still 
given by adding the electron's charge to that of a charge 
Cy = —2ps/{2ps + 1) object. In other words, the charge 
we associate with the CF, being linked to the Hall con- 
ductance and incompressibility, is unaffected by the pro- 
jection, although the notion of each electron having 2s 
vortices sitting on (or even near) it is no longer true. 
What then is this object that pairs with the electron and 
how is one to describe it theoretically ? 

The hamiltonian theory described here provides an an- 
swer. In this theory we enlarge the Hilbert space to in- 
clude additional degrees of freedom, accompanied by an 
equal number of constraints. These new degrees of free- 
dom (prevented from having any density fluctuations by 
constraints) will turn out to have charge e^ and pair with 
electrons and change e to e* . / shall refer to them as vor- 
tices, for want of a better name, but in view of what was 
said above, they are not related in any simple way to zeros 
of -9 LLL- 

We now review the hamiltonian description, which, 
not surprisingly, is relies heavily on earlier work. 

• Lopez and FradkinEII took the first major step 
and attached 2s flux tubes by the singular gauge 
transformation-of the wavefunction (due to Leinaas 
and Myrrheimca), from electrons to Chern-Simons 
fermions: 
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and opened up the field theoretical description of 
the Jain states. This was applied to, the gapless case 
V = 1/2 by Kalmeyer and_Zhangc2l, by Marston et 
al who used bosonizationEj, and in a verji-exhaus- 
tive treatment by Halperin, Lee and ReadE2l(HLR) . 
( Flux attachment had been done earlier to con- 
vert electrons to the pjmposite boson by Zhang, 
Hansson and KivelsonEiJ and earlier still to anyons 
by Fetter, Hannah and Laughlii£3. Attachment of 
vortices instead of fiux tubes by a complex gauge 
transfopnation was explored by Rajaraman and 
SondhiN.) 

Murthy and I introduced collective coordinates, to 
describe long wavelength density fiuctuations, as 
did Bohm and Pinea^J in their treatment of plas- 
mons. For every extra degree of freedom so intro- 
duced, there was a constraint on physical states, to 



keep the problem same as before. The collective 
coordinates corresponded to oscillators at the cy- 
clotron scale. Putting them in their ground states 
and projecting to the physical sector using con- 
straints, was seen to produce the zeros that turned 
flux tubes into vortices, i.e., produced the Jastrow 
factors. 

• To expose the low energy physics, Murthy and I in- 
troduced an additional unitary transformation that 
decoupled the oscillators and the fermions. This 
was however done approximately: 

1: We worked at long distances. Thus if any quan- 
tity had an expansion in powers of ql, we kept just 
the leading term. 

2. When the density operator was encountered in 
a product with other operators, we used the RPA: 
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We made the first approximation so that we could 
introduce a small parameter ql where there was 
none. The second ensured that the operators at 
small q like p{q), did not mix with those at high q 
in the unitary transformation. These were the min- 
imal assumptions we had to make before we could 
carry out the decoupling transformation. Despite 
these approximations a reasonable quantitative and 
qualitative description emerged. 

The long-distance, low-energy theory Murthy and I 
derived was given by the following set of equations: 



H = V^l f ^p(q)^;(q)p(_q) 
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= x(q) [Physical State) 
n* = p + eA* 

A* = 
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(31) 
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Note that physical states are to be annihilated 
by the constraints x- The magnetic moment of 
e/2m on eact-particle, predicted by Simon, Stern 
and HalperirH, that arises naturally here is not 
shown, and neither is the contributions to H or the 
charge from the oscillators, which are frozen in their 
ground state. The kinetic energy of the fermions is 
quenched in the small q sector if the number of 
oscillators equals the number of particles, i.e ., Q , 
the largest oscillator momentum obeys Q = \/4TTn. 



With this choice of Q, H reduces to the electro- 
static interaction between electrons written in the 
new basis. 

While this formalism is good only for small ql, it 
can still be useful. For the Zhang-Das Sarma j^DS) 
potential with A > 1 we were able to calculatea gaps 
that agreed well with the results of Park_and Jain. 
In a collaboration with Park and JainEj we also 
established some scaling relations between fractions 
at the same p (number of filled CF Landau levels) 
but different 2s (number of vortices attached) that 
seemed to work very well. 

• While the small ql theory had its share of numerical 
successes, it had some disturbing conceptual prob- 
lems. For example the constraints did not close 
to form an algebra and did not commute with the 
charge or the hamiltonian built out of it, except to 
leading order in q. This implied that charge was not 
gauge invariant and the physical sector not defined: 
how could constraints at two different g's annihi- 
late the physical states but not their commutator? 
While gauge invariance can be implemented order 
by order in a coupling constant, this is not so with 
respect to q which is integrated over. (One can use 
Q as a small parameter—but some of the central 
physics gets modified.EZl-t3) 

These pijoblems were resolved in my minimal 
extensiontl of these results to all ql that is mathe- 
matically and physically attractive. Let us assume 
that Eqns. ( |j^ - pq j represents the beginnings of two 
exponential series and adopt the following expres- 
sions for charge, constraint, and hamiltonian: 



p(q) = Y. exp(-*q • (r, - ^^z x H*,)) (38) 

3 

= Y^ e-*'i-^-J (39) 

3 

P 

X(q) = Y exp(-*q • (r, + z x n%)) (40) 

3 

= Y e-"i-^-J (41) 

3 

H=^ f J^p-(q) i;(g)e-('')^/2 p(-q) (42) 



2 J (27r)2 

Note that Re and Ri, were fully determined by the 
two terms we did derive. The gaussian in Eqn. (W3) 
will be explained shortly. 

This is my final answer. 



III. ANALYSIS OF THE CF HAMILTONIAN 

Readers who either skipped the derivation or were 
troubled by the approximations, are invited to take Eqns. 
(pSfji^) as an efi'ective field theory of CF for which the 



author can provide a plausible lineage going back to the 
electronic hamiltonian. 

To understand what these equations imply, let us be- 
gin with the coordinate appearing in the exponential in 
the expression for p{q): 



Rp = r 



Its components obey 



f 



(1 + c) 



z X n* 



[^ex -, -^ey\ — ^^ 



(43) 



(44) 



the commutation rules of the guiding center of a unit 
charge object. This, together with the fact that it enters 
the electron density, tells us it describes, in the CF basis, 
the guiding center coordinates of the electron. 
Next consider the coordinates appearing in x(g) 



R„ 



c(i + c; 



zx n* 



y^vx ■} ^vy\ — ^^ I C 



(45) 
(46) 



These describe, in the CF basis, the guiding center co- 
ordinates of a particle whose charge is — c^ = e^ = 
-{2ps)/{2ps + l). 

This is exactly charge of the 2s-fold vortex. It will be 
seen to pair with the electron and reduce the charge down 
to e* . For these reasons we shall refer to it as the vortex, 
although the nomenclature is not ideal. For one thing, we 
know that this object does not correspond to the 2s-fold 
zeros of the LLL wavefunction, which does not generally 
have 2s non-Pauli zeros per electron. In addition, the as- 
signment of physical meaning to objects that appear in 
an enlarged Hilbert space is at best schematic. Recall the 
oscillators, which when put in the ground state and pro- 
jected to the physical sector, produced the Jastrow factor 
with its vortices. We could, for this reason, call them vor- 
tices. But we must not forget that prior to projection, 
neither the oscillators, nor their wavefunction had any 
meaning in electronic language, and that conversely, in 
the electronic Hilbert space, there were no independent 
degrees of freedom corresponding to vortices, which are 
really made up of electrons. The whole idea of going to 
an enlarged space, as in the case of Bohm and Pines, is to 
be able to handle, in intermediate stages, collective vari- 
ables as canonical coordinates independent of electrons. 

Finally 



[Re , R^,] — 0. 



(47) 



Thus the four dimensional fermionic phase space has 
yielded two independent sets of canonical coordinates, 
Re and R^. Now this is twice as many coordinates per 
particle as in the electronic LLL problem. But the 
constraints Eqn. (^) tell us the density formed out of 
Ry has no fluctuations, so that the number of indepen- 
dent coordinates matehes the LLL. This is reminiscent 
of Bohm-Pines theoryEJ, where, once plasmons are intro- 
duced at small q's, the fermions are not allowed collective 
density oscillations at these g's. 



*|2 



zxn* l 

o'^v= •■ +c{1+C) 




zxnl 

(1+c) 



*|2 



FIG. 1. Anatomy of the CF: note that a CF at a point r is 
flanked by the vortex and electron. They are bound by terms 
in H that grow with D*. 



Consider Eqn. (|3| and |5|). They show that a CF at 
r with kinetic momentum 11* , is flanked by the electron 
and "vortex" within a distance of order PH*. (See Fig- 
ure |l|). Its total charge is their sum e* = l/{2ps + 1). Its 
dipole moment (in the frame r = 0) is d* = —PzxU*. Its 
size PH* ~ I near the Fermi surface, making it a well de- 
fined object in this energy range. The vortices are bound 
to the electrons, since as their separation (proportional 
to n*) grows, so does the energy since H will be seen to 
have terms that grow with 11* . Thus we have managed to 
reach one of the goals of any theory of the FQHE: explain 
how CF get their kinetic energy from the electrostatic en- 
ergy of the electrons. Indeed the entire CF hamiltonian 
is just the electrostatic interaction of electrons written in 
the CF basis. While these features are interesting, they 
are still heuristic. First, the discussions involving oper- 
ators are semiclassical. Next, 11* is not a constant of 
motion except dX v = 1/2 when it equals p. Finally the 
first quantized expressions do not include the effects an- 
tisymmetrization. However, as we proceed, we will find 
an operator expression of these ideas that is robust and 
survives in a second quantized theory fermions. 

Given the commutation relations of Re and R„ we 
easily deduce those of p and x'- 



[p(q),p(q')] =2isin 

[x(q),x(q')] =-2isin 
[X(q),p(q')] =0. 



(q X q') I 



l\ 72 



p(q + q') 



(48) 



(q X q')P) 



'^;2^^ 



2c2 



X(q + q') (49) 
(50) 



One sees that p is not only algebraically closed, but 
obeys the Girvin-MacDonald-Platzman GMP algebraEj 
for magnetic translations. I will keep referring to it as 
the projected charge density, from which it differs by a 



factor e ("^'^ /''. In the hamiltonian Eqn. (p^), which is 
just the electrostatic interaction written in the CF basis, 
the factor e^^'''^ /^ takes this difference into account. 

Note that p constitutes a specific reahzation or rep- 
resentation of the GMP algebra in terms of the final, 
canonical, CF variables, a feature that allows one to ap- 
ply standard many-body methods to H . 

Finally the constraint commutes with the projected 
charge and hence the low energy hamiltonian which is 
quadratic in it. As in the Yang-Mills case, the constraints 
form a nonabelian algebra and commute with H. 

How does the small q - RPA treatment manage to gen- 
erate these coordinates i?e and i?„ with their nice com- 
mutation relations? This is my current understanding. 
The first two terms in p and Xi which fully determine 
these coordinates, could be derived in another theory, 
equivalent to ours for small ql but not all ql: the theory 
explored by Stern et aH. In this theory Qj the upper 
cut-off for the oscillators is, assumed to be vanishingly 
small. This means that we can safely assume that ev- 
ery q in the problem, including those integrated over, are 
small, being bounded by Q. Next, RPA becomes exact 
in this case since non RPA terms involve multiple q inte- 
grals which are suppressed by higher powers of Q in the 
unitary transformation. In such a derivation one gets ex- 
actly the same first two terms. These encode the charge 
and dipole moment of the CF, which characterize the 
CF in the infrared. While my extension to all q is math- 
ematically satisfactory, it need not be numerically exact 
down to arbitrarily small length scales. We should have 
been prepared for this since one cannot have at the same 
time a scheme that is analytically tractable and numer- 
ically exact, unless we are dealing with exact solutions. 
We managed to make the nonperturbative passage from 
electrons to CF by the exponentiation, which in turn was 
guided by the LLL algebras. 

Recall that all proofs of nonzero static compressibil- 
ity at V — 1/2 relied on a care£uJ. irtinlementation of the 
constraints or gauge invariancctHJljJlj. We may now un- 
derstand this as a follows: x = means that only the elec- 
trons in the CF respond to the static potential, exhibiting 
nonzero static compressibility of unit charge objects. 



A. How to solve H? 



Now we ask how we are to handle Eqns. (|3§|-[42D. As 
shown in Appendix 2, p-filled LL's and particle and hole 
excitations on top of it, arc HF states of our H. (The 
proof relies on the rotational invariance of the hamilto- 
nian.) One possibility is to ignore the constraints alto- 
gether and proceed with the HF approximation. This 
will however lead to the following fatal flaw: transition 
matrix elements of p will start out as q and the structure 
factor S{q) will go as q^, in violation of Kohn's theorem. 
(The the q^ sum rule is saturated by the oscillators that 
were decoupled.) 

We must therefore bring in the constraints and hope 
they will bail us out. A standard way to incorporate first 



class constraints is to introduce them into the path in- 
tegral with a Lagrange multiplier and try to solve the 
theory in an approximation that respects the gauge sym- 
metry i.e., the constraints. We will discuss it shortly, but 
in the present case I will use a solution that is essentially 
what Murthy and I used in the small q theory.0: replace 
p by the preferred combination for charge density 



p^ = p-c fx 



where 



/ 



^ g-9';V8ps 



Vortex Form Factor 



(51) 



(52) 



This combination is equivalent to p in the physical sector. 
The factor e^'^ ' /*p* (absent in the earlier work and 
unimportant for any fraction other than 1/3) takes into 
account the fact that since the vortex and the electron 
have different magnetic lengths, to convert the vortex 
magnetic translation operator to the magnetic number 
density we need a different gaussian from the one we 
absorbed into H in Eqn. (|4^). Note that p^ is weakly 
gauge invariant. 



[X , P^] - X- 



(53) 



Clearly so is the H{ff) that I shall use. Weak gauge in- 
variance is enough to keep physical and unphysical states 
from mixing. 

Consider the series expansion of ff: 



r-E' 



-jqr, 



2ps 



— -ii^qx n*j+o- (qx n*jY 



(54) 



If we expand e~*''""j to first order in q we can verify 
that the term linear in q contains only the guiding center 
coordinate of the CF (r — l*'^z x 11*) with no admixture 
of the cyclotron coordinate. Thus it does not contribute 
to the order q transition matrix element. This is the 
unique multiple of the constraint we can add to p, with 
this property. 

With the constraint implemented this way, we are in 
compliance with Kohn's theorem. But there is more. 
Consider the series Eqn. (p4). The first term, propor- 
tional to CF density has the coefhcient e* . The next term 
has the dipole moment given by Read using wavefunction 
arguments. (This was done for i/ = 1/2 and is expected 
for the whole series, as suggested by Figure 1.) The van- 
ishing of the third order term explains the success of the 
small q theory. 

It can also shown that if only terms linear in q are kept 
in ff, the algebra closes with 



(55) 



sin(q X q' 1^/2) ^ q x q' P/2 



in the structure constant of the GMP algebra Eqn. (|4S 
The significance of this is not known. 

It is remarkable that a single guiding principle, Kohn's 
theorem, leads to a combination with all these properties. 
Since the internal structure of the CF is built in at tree 



level, we expect that vertex corrections (due to the con- 
straints) must vanish as q ^ 0. We shall not refer to the 
constraint any further. 

Recall the cautionary note about ascribing meaning to 
objects in an enlarged Hilbert space. Why should we give 
any significance to the terms in Eqn. (p4)? After all, we 
could have added any multiple of the constraint x to p in 
the physical sector and each would have had a different 
monopole and dipole term in the series expansion. This 
is indeed true and in an exact calculation it would not 
have mattered how much of the constraint we added, the 
final e* and d* would have been the same. However in 
our HF calculation we have included the constraint in an 
unusual way, guided by Kohn's theorem. That is, we have 
implemented the constraint to the extent we ever will. 
Since this led to a unique preferred combination (at small 
q) , we assume that any corrections due to constraints will 
affect only the higher order moments. 

In the operator approach there is no problem with how 
to assign this or that vortex to an electron. All one claims 
is that the density operator that obeys Kohn' theorem in 
our HF calculation couples to any external potential with 
a charge e* and a dipole moment d* = P qx 11*. 

There are no problems with antisymmetrization 
among particles: we simply express this first quantized 
density operator in second quantized form with Fermi 
operators. 

Let us return to the standard way of incorporating 
first class constraints in an approximation that respects 
the gauge symmetry i.e., the constraints. While this has 
not been done for general fractions, it has been done by 
Readuj for bosons at ly — 1, which turn into fcrmions in 
zero field uppa single vortex attachment a la Pasquicr 
and Haldanetd. (The lessons learnt from this exercise 
are directly applicable to us since the only difference is 
in the coupling constant of the gauge and matter fields.) 
Implementing the constraint in a conserving approxima- 
tion leads to a propagating gauge field whose longitudinal 
part ql screens the charge fully, leaving behind dipoles of 
moment d* which then interact via the transverse gauge 
field ax- The propagator of ax, just as in HLR, favors the 
region iuj ~ q'^ . (The field a^, also produces mass diver- 
gences at the Fermi energy as in HLR and restores com- 
pressibility.) However, away from the ultra-low frequency 
region, the answer is given by the correlation function of 
independent objects of charge e* = and dipole moment 
d* = P qx p. 

For gapped fractions (not too close to i^ = 1/2) and/or 
at T > 0, a description in terms of independent particles 
with the right e* and d* is likewise expected to be a good 
approximation, likely for all co, since the gap and/or T 
will cut-off the low frequency end where ot raises its 
head, and the major effects of a^ are already encoded in 
e* and d*. Rather than reach this description using the 
conserving approximation (which is pWry difficult away 
from V — 1/2 due to LL structure,c3) I use a scheme 
Murthy and I devised for small ql. (Away from v = 1/2, 
ot could affect the statistics of the quasiparticles.) 

For the benefit of the readers who just joined in, I dis- 
play the equations will be used in the subsequent calcula- 



tions: 



H' = \j ^r (q) «(g)e-('')'/2 r (-q) (56) 



r (q) = P(q) - c'e 



2-qH^/»ps 



x(q) 



p(q) =^^^^{-^^-{^j-Yi;~J' ^ n*j)) 



(57) 
(58) 



X(q) =^exp(-iq-(rj-+ _^ z x n%)) (59) 



[H^,X] ^X 

— x\ physical state) 



(60) 
(61) 



Many leading long wavelength effects of the constraints 
have been built into ff; they will be ignored in the subse- 
quent HF calculations, as will be the superscript on H^ 
since we shall always use this expression in terms of the 
preferred charge ff. 

We separate H into free and interacting parts 



H = Ho+Hi 



(62) 



wherein the two pieces Hq and Hj correspond to diagonal 
(i = j) and off-diagonal {i ^ j) terms when the double 
sum over particles is expanded. In the simplest case v — 
1/2 we have 



^o(^ = l/2) = E2/S 



sm 



q X k,P 



v{q)e' 



(63) 



If we expand the sin in a series and keep the lowest term, 
we get an expression quadratic in momentum from which 
we can define an effective mass 1/m*. Thus we have man- 
aged to generate the CF kinetic energy,m operator form, 
in terms of the electron-electron interaction. But we have 
more. First, there are more powers of momentum in the 
kinetic energy. (We shall however see that only quartic 
term is important in HF.) Next, there is also Hj, which 
can modify all the numbers. Thus l/m* depends on the 
momentum and will be defined at tbe,Fermi surface. For 
J/ = 1 bosons, Haldane and FasquierEll obtained the same 
Hq by algebraic methods aimed at a direct LLL formal- 
ism. In the present case we have H for the entire Jain 
sequence. Away from ^ or i, H is even more compli- 
cated, but of a very definite, known, functional form. 



IV. THE ACTIVATION GAPS OF FULLY 
POLARIZED STATES 

Here we use the hamiltonian theory to compute acti- 
vation gaps for fractions i^ = p/{2ps -I- 1) in a field so 
strong that the system is fully polarized. We will probe 
the theory in the following ways: 

• We will compare the gaps to ihose obtained by 
Park, Meskini and Jain (PMJ)I3 using trial wave- 
functions for the Zhang-Das Sarma (ZDS) poten- 
tial 



v{q) 



27re2 



-ql\ 



(64) 



for p — 1,2,3 and 4 and s = 1. These serve as a 
benchmark, at least for fractions not too close to 
i/= 1/2. 

• For the benefit of other users, a fit to the gaps in 
the experimentally significant region 1 < A < 2 will 
be given for both s = 1 and s = 2. The gaps are 
also expressed in terms of an effective activation 
mass 771a. 

• We willjCompare the theory to the experiments of 
Du et aM and Pan et aH. 

• We will examine the charge density profiles in states 
with a quasiparticle or quasihole. 

• The gaps are computed for a gaussian potential and 
compared to PMJ for p = 1, 2, 3, 4 and s = 1. 

• The effect of turning off the interaction Hi in H = 
Hq + Hj, will be explored. 

• We will compare the resulta,to the exact diagonal- 
ization results of Morf et aB. 

• It will be seen how well particle-hole (PH) symme- 
try works. For example, is the activation gap for 
2/5 same as that for 3/5 in the fully polarized case? 



A. Comparison of gaps in the HF approximation to 
PMJ 

We use the HF approximation. As shown in Appendix 
2, p-fiUed CF LL and particle-hole excitations thereof are 
HF states of our H. The HF ground state is given by CF 
filling p LL's, which we will denote by |p). I will use a 
boldface symbol such as p to label a Slater determinant 
with p occupied Landau Levels. Nonboldface symbols 
will label single particle states. Note also that the actual 
LL index n for the state labeled hy p is n = p — 1 since 
the LLL has index n = 0. 

The gap is defined by 



A = (p + PH\H\p + PH) - {p\H\p) 



(65) 



where PH stands for a widely separated particle-hole 
pair. Now that this is exactly hpw gaps are computed 
in PMJ's wavefunction approachQ. There is however one 
big difference hidden in the notations. There too the 
hamiltonian is just the interaction, but written in the 
electronic basis (with p{q) — ^ ■ exp(iq • r ■)) while the 
states (which carry the same label) are these simple wave- 
function, multiplied by the Jastrow factor and then pro- 
jected to the LLL. Projection leads to a very compli- 
cated expression for the wavefunctions. In the present 
approach we have tried to incorporate these effects by 
going in the reverse direction, from electrons to CF's, 
and obtaining complicated expressions for the charge and 



other operators, but with simple expressions for the wave- 
functions. While these operator expressions are unusual 
in form, they are still simple enough in appearance and 
amenable to exact analytical treatment, because of the 
approximations that went into the derivation. Thus does 
not expect the present results to match those of the Jain 
approach in their accuracy. This is indeed the case, un- 
less the potential is fairly soft (A is larger than, say unity). 
The aim of the present approach is to provide a 10 — 20% 
theory for soft potentials (which do seem relevant to ex- 
periment) in which the calculations can be performed 
analytically and the physics of the quasiparticles is trans- 
parent. 

Rather than work with a widely separated particle- 
hole (PH) pair, I first find the energy in a state with just 
the particle and add to it the energy of a state with just 
the hole and subtract double the ground state energy. 
While the details are relegated to Appendix 4, here is 
the central idea. 

One begins with the second quantized expression for 
the preferred charge operator p(q): 



'(q) = E 4 



m2n2 ^mini Pm2n-2 ;mini 



(66) 



ni2n2:mi7ii 



where d\^^^^ creates a particle in the state \m2 n2) where 
m is the angular momentum and n is the LL index of CF 
in the weakened field A* — A/{2ps + 1) with a magnetic 
length 



I* = l^/2ps + l. 



(67) 



The key ingredient in the HF calculation is the matrix 
element Pm2n2;mini- It is shown in Appendix 1 that this 
matrix element factorizes: 



Pm2n2\mini — Pm2mi **" Pn2ni 



(68) 



The gaps depend only on p" „ , the superscript on 
which will be generally dropped. 
As shown in Appendix 4, 

/j2 
^ exp(-a;/(2ps + l))viq)G{p, x, s) (69) 

_q^l*^ _q^P{2ps + l) ^^^^ 

(71) 

(72) 
(73) 



2 2 

G{p, X, s) = Go{p, X, s) + $(p, x) 

Goip, X, s) = c^e""^^/^ exp ' 



\{Aps){2ps + I) , 

X {Lp^l{xT) - Lp(xT)) 

p-i 
*(P: x) ^^{R (P- 1' Tn,x)- R {p, m, x)) 



i?(a,6,x)=x"-''(- 



-c^-'-^exp 



Lt'ix/c^) 



( - 

\Aps{2ps+ 1) 



2ps 



2ps + l 

T ^{C~ 1/c) 



(74) 

(75) 
(76) 
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• At A = 0, the coulomb case, the gaps are finite 
in contrast to the small q theoryQ. This is due to 
the gaussian factor e^'^ ' ^^ which was absent there. 
The slope of the graphs in the present theory is 
nonzero at this point. It is readily verified that 
dA/dX at A = is the gap due to a delta- function 
potential, and should vanish for spinless fermions. 
It does happen for PMJ, whose electronic wave- 
function is explicitly antisymmetric. The present 
theory of CF does not give good answers for poten- 
tials as short ranged as the delta function. Indeed 
even the coulomb interaction is too singular, and 
the theory begins to work well only beyond A ~ 1. 

• Beyond A ~ 1 the agreement is quite fair in general 
and very good for 2/5. 

• The gaps do not vanish for any fraction and any 
finite A. 



0.0 0.5 



1.5 2.0 2.5 3.0 



FIG. 2. Comparison of dimensionless activation gaps Sa to 
the work of Park, Meskini and Jain ei al for the fractions 
1/3, 2/5, 3/7, 4/9 (p = 1, 2, 3, 4, and s=l) as a function of A, 
the thickness parameter in the ZDS potential. 



where L^ is the associated Laguerre polynomial. 

Often we will use the dimensionless gap 5 defined by 



el 



It is useful to know that 

Aa e2 Sa 



ks elks -^«v/^^"°^ 
I shall use the values given in Table (0): 



(77) 



(78) 



eB 
rriekB 

elkg 

e 


1.34 B{TYK 


50 ^BiTYK 


.026 me^B{T) 



Table 1 Approximate numbers used in this paper, with 
ks the BoUzmann's constant, B(T) the field in Tesla. 



Figure g shows the gaps computed for 1/3, 2/5, 3/7 
and 4/SI.for the ZDS potential and compared to the work 
of PMJ0 in the region < A < 3. The following features 
are noteworthy. 



B. Activation masses 

I have computed gaps for many other fractions, in- 
cluding for 3 = 2, when four vortices are attached to 
form CF's. 

Rather than show more plots, I will now analyze the 
theory in terms of nia the activation mass defined by 



A„ 



eB* 
rua 



eB 



{2ps + l)ma 

Comparison to Eqn. { \n\ j shows that 

1 e'^l e^l 

— = — <5,(2ps + l)^— C,. 
ma e e 



Thus 



Ca=5a{2ps+l). 



(79) 



(80) 



(81) 



The significance of Ca is that it approaches a limit 
as we_approach v = 1/2 or 1/4, as first emphasized by 
HLREil. What we will see now is the Ca does indeed have 
a nice limit, but this limit depends on A, a parameter 
that was set equal to zero (coulomb case) in HLR. We 
will focus on the s = 1 series, which converges to i^ = 1/2 
as p -^ oo. To remind us of this fact a superscript 2 = 2s 
will be appended. Consider the Table (0), of gaps fitted 
as a function of A in the interval 1 < A < 2, although the 
fit will work for modest excursions on either side. 



p 


Ai'Vfcs = 50v/s(r)5r 


S^^^ 




1 

2 
3 

4 


5.31v/S(T)/A 
2ml/B{T)/X 
1.23^B{T)/X 
0.87y/B{T)/X 


.106/A 
.042/A 
.025/A 
.017/A 


.32/A 
.21/A 
.17/A 
.16/A 



Table 2 Activation gaps as a function of X for 1 < A < 2 

according to the hamiltonian theory. Note the conver- 
gence of Ca as p ^ oo. 
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Table 3 Comparison of activation masses to Du et al, 



Based on Table (^, and a similar one for states near 
V ~ 1/4, I find that we may write, near these states, 






.160 
.148 

A574 



(82) 
(83) 



where the expression, including the exponents (1, |) are 
approximate. I— . 

Consider the normalized mass defined by Pan et aO 



m„ 



.VW) 



(84) 



where m^ is the electron mass and B{T) is the field in 
Tesla. 

In terms of C, 



.026 



Combining Eqns. (p2[pq) 



.163A (s = 1) 
.175A5/4 (.s = 2) 



(85) 



(86) 
(87) 



We find that the suitably scaled masses rrinor are com- 
parable for s = 1 and s = 2 but not exactly equal. No 
fundamental arguments exist for their equality since the 
answer depends on the potential, parameterized by A. 



C. Comparison to data 

In comparing these gaps to experiments, I will limit 
myself to z/ < 1/2. States like 3/4 are related by particle- 
hole symmetry if full polarization is assumed, and states 
with v > I require assumptions about filled electronic 
LL's which I do not want to make. ■— ■ 

Consider the experiments of Du et alllJ, who have ex- 
tensive data on activation gaps. Given that the experi- 
ments, unlike PMJ, have an unknown contribution from 
LL mixing and impurities, it is not clear how to apply 
the theory. There is no ab initio calculation that in- 
cludes these effects. (There is however reliable evidence 
that LL mixing is arvery small effect at the values of A 
under consideration.E^). I will therefore compute gaps 
using the ZDS potential with A as a free parameter, and 
ask what A fits the data, just to get a feel for its size. 
The results are summarized in Table y. 



v 


B{T)) 


Ar^cx) 


Af™(°/s:) 


A 


1/3 

2/5 
3/7 


13.9 
11.6 
10.8 


8.2 
3 
2 


5.3VS(r)/A 
2.08v/S(r)/A 
1.23y^B{T)/X 


2.4 
2.4 
2.0 



sample A, which has a density n = 
last column gives the best fit to A. 



1.12-10''cm- 



The 



Is it possible to describe the disordered sample by 
some effective A? It does not seem likely, given these 
three data points: the line through them (plotted against 
B* = B — 9.27 Tesla) has a negative intercept, while the 
pure system calculations give a gap that never vanishes 
for any finite A. Furthermore, the negative intercept is 
between 1° — 2"K, while the gaps are at best 6 — 8°K. 
Thus the effects of disorder appear to be quite signifi- 
cant. Therefore no attempt will be made to find an effec- 
tive A. If one day we get samples for which the disorder 
broadening is much smaller, we can attempt this. In the 
meantime, LDA and exact diagonalization calculations 
suggest that the answers for the pure. Sjystem differs by 
roughly a factor of two from the datal3,El. 

Consider now the results of ParJid et al. Rather than 
attempt to fit their gaps to the theory, I consider the 
following issue they raise about the normalized activation 



m,^B{T) 

They observe that these masses arc in the range .25 — .35 
near v — 1/2 and 1/4. How does this rough equality of 
masses of fermions with two and four vortices fit in the 
present theory? If we compare their results to Eqns. ( pq - 
^) we extract the range of values for A listed in Table 
r They seem to be spread over a range that is more or 
less equal. The theoretical prediction for this case is sen- 
sitive to how exactly the two fractions are reached. For 
the present case, wherein n is fixed and B is altered, we 
expect A'^'^V-^*'^'' — V2- I suspect that the experimental 
A's do not show this factor of V2 because the effects of 
disorder are most pronounced near the gapless states. 



Theoretical value 


A implied by data 


rn^"" = .163A(2) (s = 1) 
m^'^ = .175(A(4))5/4 (s = 2) 


A(2)=1.5-2 
A(4) = 1.3-1.75 



Table 4 Given that normalized masses for 2 and 4 vor- 
tex attachment are in the range .25 — .35 the table asks 
what the corresponding values of A are. Theory predicts 
that A(4) = V2A(2)s. 



D. A study of particle and hole profiles. 

Now we consider the charge densities in a state with 
either a widely separated particle-hole pair in one of the 
gapped fractions or just a CF aX v = 1/2. The calcula- 
tions are detailed in Appendix 3. 

Figure (0) which shows the distribution of charge in 
units of the ambient density, for three gapped fractions 
with a widely separated particle and hole. To compare 



12 



L.3 


1/3 






f 




-\- 






II 


/ 5 


10 


15 


-nrV 21 


3.9 
3.8 
3.7 


/ 










approximation where the function drops to about half the 
ambient density. This breakdown of the present picture 
at short distances is generic. Note also that while our 
unitary transformations attempt to go from electrons to 
CF while Jain's wavefunction approach goes from CF to 
the electrons, the two are not of equal accuracy. Given 
the simplicity and tractability of the final operators in 
the present model (which cannot possibly carry the com- 
plexity of the projected Jain wavefunctions in electron 
coordinates) this is to be expected. 



E. A closer look at 



1/2 



Kt V = 1/2 the CF of mopaentum k is predicted to 
have a dipole moment I'^z x knJ. Let us understand this 
statement. It stems from examining the trial state 



* 



^n(^ 



ZjfDet\e 



(89) 



i<j 



Before projection, the vortices are on the electrons. The 
CF charge is e* but the dipole moment is zero. To 
project, one begins with 



Particle and hole profiles 



ikr 



exp —{kz + kz) 



(90) 



FIG. 3. Particle-hole profiles, in units of the ambient den- 
sity, at 1/3, 2/5, 3/7. The results are compared to the unpub- 
lished work of Park and Jain by placing them at antipodal 
points of a sphere, with the hole at the north pole r = 0. The 
coordinate r is the distance on the sphere measured in units of 
I. The solid lines are my predictions, the horizontal (vertical) 
lines show the heights (locations) of the maxima and minima 
observed by Park and Jain. 



to the unpublished work of Park and Jain who placed 
them at antipodal points of a sphere, I have done the 
same. The hole is placed at the north pole r = and 
r is the distance along a great circle in units of I. The 
solid lines refer to my results. The small horizontal lines 
indicate the values of the maxima and minima of Park 
and Jain results, while the small vertical lines show their 
locations. Three features are noteworthy. 

First, the density is best predicted at 2/5. It is not 
surprising that the gap follows PMJ best in this case. 

Next, in the case of 3/7, the particle and hole actually 
overlap and get entangled in the Monte Carlo work on 
a finite sphere. In the present infinite volume approach, 
we have no trouble distinguishing them apart since they 
were computed individually and superposed to draw the 
figure. 

Finally, consider the case of 1/3, where one expects 
to find the best results for gaps (due to a large gap), 
but does not. Since there is just one filled CF level, if 
you make a hole at r = the density must vanish. This 
feature found in Park and Jain, is absent in the present 



(where k = k^ + iky, ) recalls that z -^ 2Pd/dz, and 
sees that the derivatives in the exponential will move z 
to z -|- ik. This suggests that the vortex will move off 
the particle by an amount ikP, which explains the origin 
of the dipole and its moment. This analysis in terms of 
just one particle corresponds to what was inferred from 
Figure 1 or the analysis of Re and R„ . 

While this is true, it does mean that the spatial distri- 
bution of zeros in the LLL wavefunction (as a function 
of one coordinate, the others being frozen) will resemble 
the momenta that fill the Fermi sea. There are several 
reasons for this. 

• Besides the one exponential we considered, there 
are other exponentials that move every other coor- 
dinate by the corresponding fc, so that in the end 
{zi — ZjY —^ {zi + iPki — Zj — iPkjY. This de- 
pendence on the difference of momenta is in ac- 
cord with the symmetry first noted by-nHaldane, 
and dubbed K-invariance by Stern et a£3. 

• All of the above refers to just one term in the ex- 
pansion of the determinant. We must still antisym- 
metrize over all possible pairing of momenta with 
coordinates. 

• Many of these zeros will have to disappear upon 
projection, leaving just a total of two per particle. 
It is far from clear how the zeros will be distributed 
in the end. All we know is that one has to lie on the 
electron by the Pauli principle, and thejOther forms 
a single vortex, as in D.H. Lee's picturea. The lat- 
ter must be at a distance 2kP from the electron 
to preserve the dipole moment. Since the "size" of 
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the CF at the Fermi energy, 21, is close to the inter- 
particle spacing of 2-^7?/ ~ 3.5/, there is definitely 
going to be some ambiguity in pairing zeros with 
electrons. 

These complications notwithstanding, the dipole 
makes a clear appearance in the operator approach. It is 
not tied to zeros of the wavefunctions and emerges as fol- 
lows: if the expression for the preferred charge density ff , 
(which obeys Kohn's theorem), is coupled to an external 
potential, the second term in the series Eqn. (m) is pre- 
cisely that of a dipole of strength Pz x k. This operator 
definition is unchanged by antisymnietrization (we sim- 
ply write Eqn. ( |54| ) in terms of fermion operators in sec- 
ond quantization) and crystallizes the notion of the dipole 
inspired by an analysis of the trial wavefunction. The 
dipoles also emerge in the high frequency density-density 
response computed in the conserving approximatianE3, 
as stated earlier. The zeros of trial wavefunctiongc3 do 
not seem to be the likely place to look for evidence of 
dipoles. 

I will now illustrate these points in a many-body state 
with an extra CF. It however takes a bit of work. Suppose 
we create a CF in a state of definite momentum p above 
the Fermi surface. If we evaluate f^ir) in this state, we 
get zero. This is because the CF has an equal amplitude 
to be at all places and the dipolar density gets smeared 
out completely. If we localize it at the origin, we use 
all momenta, and hence all values of the dipole moment 
which again averages to zero charge density. Consider 
however the following state 



density 




FIG. 4. The charge density of a CF showing its dipole mo- 
ment. The CF is localized in x but spread out uniformly in 
y. Luckily this does not smear out the dipole moment. 
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FIG. 5. Comparison to the work of Park, Meskini and Jain 
et al for the gaussian potential v{q) — 2-Ke^l e"' ' " for 
p=l,2,3,4 



bo) 



dpxd\px,Po)\FS) 



(91) 



where \FS) is the Fermi sea, and po > pf- This is a 
superposition of states of all values oi Px, localizing the 
CF at X = but with fixed py = po- (As po > pp, all 
these states lie on a line outside the FS.) The particle is 
thus spread out completely in the y-direction. However, 
its dipole moment is fixed at Ppo along the a;-axis. A 
calculation done in Appendix 3 yields the following result 
in compliance with these expectations: 



p(r) 



exp 



-{x-^oiy 



exp 



'{x + \pol 



2\2 



(92) 



depicted in Figure (Q). 

If one evaluates the dipole moment of this charge dis- 
tribution, one finds it equals Pol^- This is the sense in 
which the dipole moment appears in the the operator 
approach. 

I conclude with one significant difference between e* 
and d* . The CF charge e* is robust under projection, 
while d* is not. In the unprojected wavefunction at 
V = 1/2, the dipole moment of the CF is zero since the 
vortices are on top of the electrons. The charge is zero as 
well. Consider a system with an interaction comparable 
to the cyclotron gap, such that the best wavefunction is 



the above unprojected one. Based on this wavefunction 
we would assign to the CF a charge e*, but no dipole 
moment. Consider now a problem where the projected 
state is the best. Instead of applying V in one shot, 
imagine slowly reducing the non-LLL component of the 
wavefunction to zero. Along the way, the zeros get ripped 
off the electrons and some of them disappear. Through 
it all e* is invariant, since it depends on the Hall conduc- 
tance, which is constant. But the dipole moment, which 
describes the internal structure of the CF, changes from 
zero to some value in the LLL. This value seems to be 
Z^ z X k, in the wave function based arguments of Read, 
in the operator series and the density-density response. 
I do not know if the dipole moment is robust within the 
LLL, but suspect it is, since Kohn's theorem (which lim- 
its the matrix element in the LLL) gives a unique answer 
in the operator approach. 

Some caution must be used in looking for this mo- 
ment. Even if we know how the CF couples to an ex- 
ternal potential, the response function will depend on 
the Hamiitonian as well. In particular, we know that at 
w = 0, g ^ 0, the compressibility will not vanish as q^ 
but as a constant, because of the special symmetries of 
H , a point madQ-hv Halpt;rin and Stern and discussed at 
length in Rcfs.EZl,E3 anda. Only at high frequencies will 
the dipoles behave classically. 
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to A, the full HF gap, for p = 1, 2, 3, 4 and s = 1, 2. 
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F. Other potentials 

Figure ^ shows a comparison to the PMJ results for a 
gaussian potential 



v{q) = 27re2/ e"'?'''/^ 



(93) 



Note that except for v ~ Xji the agreement is excep- 
tional. This is the kind of potential for which the present 
theory works best. 

On the other hand for a potential, 



v{r) 



(94) 



the agreement is worse than for the coulomb case since 
this potential is just as bad as r ^ and does not give 
the large r values a chance. Likewise 1/r^ fares worse 
than l/r. 

From playing with these potentials and using the PMJ 
results as a benchmark we can thus learn when the 
present model can be trusted. 



G. Effects of CF interactions 

Figure shows what happens to the gaps if the CF 
interactions are turned off. These correspond to the con- 
tribution from the Gq term in Eqn. ([7l|). Note that 
interactions seem less important for v — 1/4 and system- 
atically get less important as p increases. 

There is some freedom in defining the measure of in- 
teractions, which I took to mean the effect of Hj on 
Aa- But the effect of Hj is two fold: it renormalizes 
the self-energy of the particles and also mediates inter- 
actions between them. One can envisage a situation in 
which individual energies get strongly renormalized by 
Hi but the dressed particles are barely interacting, i.e., 
their energy barely varies with separation. In this case 
one could argue that CF are weakly interacting. Such a 
separation, between the quasiparticle and quasihole, was 
set to infinity in our gap calculations. This energy, as a 
function of separation_Js_contained in the g-dependence 
of the magnetoexcitonE3,c3 spectrum. Since the variation 
is typically comparable to the gap, CF seem to be quite 
strongly interacting by this measure as well. 



H. Comparison to the work of Morf et al 

Morf et aB have calculated the activation gaps for 1/3, 
2/5 and 3/7 by exact diagonalization of finite systems, 
paying great care in extrapolating to the thermodynamic 
limit. The potential they use is 



v{q) = 



2ne' 



Mb?) 



Erfc (qlb) 



(95) 



where b is the analog of A. Our numbers are compared 
in Figure M. 

The following features are worthy of note. 

• The calculated gaps always lie above the exact di- 
agonalization results for the two fractions shown 
(as well as for the 1/3 case, not shown). This re- 
sult agrees with the general belief that HF always 
overestimates the gaps by neglecting fluctuations. 
Compare this to the case of PMJ where for 3/7 and 
4/9 the calculated gaps were sometimes lower and 
sometimes higher. 

The fact that the theory predicts gaps that always 
exceed the exact diagonalization results suggests 
the possibility that the problem may lie not in H^ 
but in the HF approximation. If H^ were solved by 
a more accurate method than HF, the agreement 
might have been better. This however merely re- 
mains a possibility till someone solves it by, say, 
exact diagonalization. This will be hard since the 
theory has constraints and is formulated in the full 
fermionic Hilbert space. 

• The general agreement is worse for this potential 
than for the ZDS case. This is because at large q 
this potential goes as 1/q while the ZDS potential 
falls exponentially. 



I. Particle-hole symmetry 

For the fully polarized case one expects that within 
the LLL, Aq for v and 1 — v will be equal. Let us focus 
on just the case s = 1, which corresponds to two- vortex 
attachment. The fractions related are now v and 
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l-I/ = 1- 



p 



p + 1 



2p+l 2(p + 1) - 1 



(96) 



Thus to find the gap at 1 
1/ = p/i'2p + 1) we must 



v corresponding to a certain 



• Replace phy p' = p+1 and consider fractions of the 
form p' /{2p' — 1). The flux quantum per CF is now 
— 1/p' which means the double- vortices have over- 
turned the applied field and changed its direction 
and p' levels are filled in this weakened and reversed 
field. In Eqns. ( pOffZ^ ) we must make the obvious 
changes, p ^ p' = p+l and {2p+l) ~> (2p' — 1), in 
all the denominators, and the not so obvious change 
exp(-a;/(4p(2p-H 1)) -> exp{x/{4p'{2p' - 1)) in the 
double-vortex form factor which reflects the fact 
that the vortex charge is bigger than that of the 
electron. 

• If we are working at A = ?A 7^ 0, we have to keep A 
constant in the comparison. Thus we need to verify 
that 



A(i/, A) = A 1 - I/, A 



(97) 



or for the dimcnsionless gap Sa — Aa/{e^/el)), 



Sa{i^,X) = 



l-v 



5a{l-V,X 



l-v 



) (98) 



My results are summarized in Table 0. 





6a{v.Q) 


<5a(l-l^,0) 


5a{l^.l) 


-J-Ii-'^Vt^) 




^{l-y)/u 


^(l-y)hy 


1/3 


.213 


.203 


.106 


.101 


2/5 


.097 


.098 


.042 


.043 


3/7 


.063 


.065 


.025 


.026 


4/9 


.047 


.048 


.017 


.018 


5/11 


.038 


.039 


.013 


.014 



Table 5 Particle-hole symmetry in the LLL for the po- 
larized cases with v = p/{2p-\- 1) requires that the dimcn- 
sionless gap 6a{i^,X) — \/-^—^da[l — i^,X,/-r^]. The 

Table considers A = and 1 and indicates that the theo- 
retical numbers obey this symmetry. 



Note that particle-hole symmetry holds very well even 
at A = when the absolute value of the gap is not close 
to the benchmark value set by PMJ. If one looks at the 
expressions for the gap, one sees that this agreement is 
very nontrivial since the two problems are very different 
in the CF formalism: they have different number of filled 
LL, and the matrix elements involved are quite different. 



V. MAGNETIC TRANSITIONS AT T = 0. 

Now we turn to the behavior of the spin of the system, 
which was assumed to be frozen along the applied field. 



This topic has explored within the CF approach rather 
extensively by Park and Jaintj. 

The coupling of electron spin to the applied field is 
given a Zeeman term 



H, 



2771. 



-B 



(99) 



where g = .44, rrif, is the electron mass in free space, S is 
given by 



S = nP 



(100) 



where n is the density and P is the polarization, to which 
each electron contributes ±1. 



A. Magnetic transitions in gapped fractions 

If Hz dominates, we expect the system to be fully 
polarized {P = 1). As we lower Hz, we may expect P 
to drop. In CF theory for gapped fractions there is a 
discrete set of allowed values of P. At i/ = p/{2ps -1-1), 
these correspond to states of the form |p — r, r). These 
stand for many-body states in which p—r LL are occupied 
by up spins and r LL by down spins. In Jain's approach, 
the actual wavefunction will be such a state times the 
Jastrow factor, followed by projection to the LLL. In the 
present approach, |p — r,r) is literally the state, but the 
operators for charge and spin are obtained by canonical 
transformations. For the interested reader I mention that 
flux attachment and canonical transformations are same 
for both spins. 

Since the uniform external field couples to the g = 
component of the spin density which is unaffected by the 
canonical transformations, Hz will have the same form 
in the final CF representation. 

It is important to note that even though the states 
are labeled by free particles (in a weakened field A*), the 
problem is not really free: in Jain's version the free states 
turn to highly correlated wavefunctions for electrons, and 
our case, the states may look free, but H is not. 

In any event, the allowed values of polarization are 
given by 



P = 



p — 2r 
P 



(101) 



Thus for example, when p — A, the allowed values are 
P = 1, .5, and corresponding to |4,0) ,|3,1) and|2,2). 
Our goal is to calculate the critical fields at which the 
system will jump from one value of r to the next as Hz 
is varied. Let 



E{p - r, r) = (p - r, r|iJ|p - r, r) 



(102) 



where H does not contain the energy due to Hz- This 
will be case for the single-particle and ground state en- 
ergies, with one exception which will be clearly pointed 
out. Since Hz is diagonal in the HF states which have 
definite spin, its effects can be trivially incorporated. 
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The HF calculation oi E{p — r, r), detailed in the Ap- 
pendix 5, gives 



E{p ^ r,r) ~ 

"p— r — 1 



n 

r-1 



p— r— 1 p— r— 1 



ii— ni — 712— ( 

r— 1 r— 1 



ni=0 



m— 712—0 



where p„i„2 was introduced earlier and discussed in Ap- 
pendix 1 and 



Ap{<i)p{-q)W) = X \p{q)r. 



(103) 



The critical field B'^ for the transition from r to r + 1 
is given by : 

e_B= n 
E{p-r,r)-E{p-r-l,r+l) = g (104) 

2?7le p 

where the right hand side denotes the Zeeman cost of 
flipping the n/p spins in the LL that switched its spin. 
This discussion assumes that B is perpendicular to the 
sample. If there is a tilt 9, we write 

eB'^ n 

E{p~r,r)~ E{p-r-l,r + l) = g ^- (105) 

2me cos 9 p 

When these i?^'s were calculated, I noticed the same 
remarkable regularity first noted by Park and Jainllj, 
namely that they could be fit by a theory of free fermions 
of mass TTip (the polarization mass) that occupy LL with 
a gap Ap = eB* /mp. In this case we would have 

111 7) — 2/^ — 1) 

E{p-r,r)-E{p-r-l,r+l) = ^^ '-Ap (106) 



since {n/p) spin- up fermions of energy (p — r — 1 + 5 )^p 

drop to the spin-down level with energy {r -\- |)Ap. 

Suppose we evaluate the left-hand-side of Eqn. ( |106| ) 
in the HF approximation to H and define 



- E{p- r -l,r + 1) 



Apirfef ^pEip^r^rl^,^_^^^^ (107) 
n p — 2r — 1 



Given that H is not free, there is no reason why 
Ap{r)'^'^f should be r-independent. But it is very nearly 
so. For example at p = 6, A = 1, 

g2 

Ap(0,l,2)'''=^ = —(0.00660,0.00649,0.00641) (108) 

which describe |6,0) —^ |5,1), |5, 1) -^ |4, 2), and 
|4, 2) -^ |3,3). This r-independence of the gaps was 
true for every fraction and every value of A I looked at. 
Yet I knew that H was definitely not free since the acti- 
vation gap Aq to make a widely separated particle-hole 
pair differs from Ap by factors like 2 or 4 (depending on 



A) and turning off Hj makes a substantial difference, as 
demonstrated earlier. 

I will place this result in perspective shortly, after not- 
ing that it has a counterpart in the gapless case as well. 
But first I present the HF results for (the r-independent) 
nip and Ap = eB* /mp. At and near v — \ and i, for 
.75 < A < 2, nip may be approximated by 



1 



--c^'HX) 



mi, 



eH 



,(4) 



c^Hx) 






.087 

AV4 

.120 

AVJ 



(109) 
(110) 



For fractions like 2/5, not too close to 1/2, I will use the 
actual TTip in comparing to experiment. 

The transition |p — r, r) -^ ||p — r — 1, r + 1) occurs 
when 



S^ 



2me cos 9 



(p - 2r - l)Ap. 



(Ill) 



B. Magnetic transitions of gapless fractions 

Let us now turn to the gapless fractions 1/2 and 1/4. 
The discrete labels p — r and r of the HF states that 
count the spin- up and down LL's are now replaced by 
continuous variables k±p which label the Fermi momenta 
of the spin-up and down seas. These momenta are such 
that the total number of particles equals n: 



ki 



k_p = kp = Ann 



(112) 



where kp denotes the Fermi momentum of a fully polar- 
ized sea. 

In the gapped case there were several critical fields B'^, 
each corresponding to one more CF -LL flipping its spin, 
each describing one more jump in the allowed values of 
P. In the gapless case the situation is different. For 
very large Zeeman energy, the sea will be fully polarized. 
It will not be worth including even one fermion of the 
opposite spin since the Zeeman energy cost alone will 
exceed the Fermi energy of the polarized sea. As we lower 
the Zeeman term, we will reach a critical field at which it 
will be worth introducing one fermion of the other spin 
with zero kinetic energy. At this point the energy of a 
particle on top of the spin-up sea obeys 



£+{k+p) = g 



B 



± 



2top cos 6 



(113) 



If we lower the Zeeman term further, the polarization 
will fall continuously and be determined by £±{k±p), the 
energies of the particles on top of these two seas according 
to 



£+{k+p) -£^(k^p) = g 



Bi_ 



2to„ cos 9 



(114) 



since this equation states that the system is indifferent 
to the transfer of a particle from one sea to another, i.e., 
has minimized its energy with respect to polarization. 
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Using the fact that at these fraction when the mean 
magnetic field vanishes, we deal with a very simple ex- 
pression: 



p"(q) 



d^k , „ , . ,q X k ^^, ,t 



(115) 



it is easy to do a HF calculation and obtain 



£±{k)) = 
2/0.(,)sin^ 



kx qf 



<Pk' 



n^(|fc'|)v(|k-k'|)sin2 



k' X k ^2 



where the Zeeman energy is not included, n^ is the (step) 
Fermi function for the two species, £q and £/ represent 
single particle energy (due to what was called Hq earlier) 
and the energy of interaction of this particle at the Fermi 
surface with those inside the sea, and 



«(fc)=«(fc)e-'='''/2. 



(116) 



When this result is used to compute £+{k+p) — 
£-{k-F), I find once again that the numbers fit a free 
theory in the following sense. Imagine that CF were free 
and had a mass to„. We would then have 



£+{k+F) - £-ik-F) 



i-2 _ i.2 
+F '^ — F 

2m„ 



(117) 



What I find is that the HF number for £{k+F) — £{k-F) 
may be fit to the above form with an rup that is essentially 
constant as we vary k±F *-e., the relative sizes of the up 
and down seas (which is analogous to an rUp that does 
not depend on the index r in the gapped case) and that 
this mp matches smoothly with that defined for the nearby 
gapped fraction. 

This result is surprising because we know the CF are 
not free from a variety of reasons. Indeed the HF energies 
do not have a quadratic dispersion: for example &i v — ^ 
and A = 1 



£{k 



±F) 



(eVeO 



(k±F 



fk±i 



V kF 
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where a == .075, 6 = -.030. 

The proper interpretation of this free-field behaviour 
will be taken up shortly. For now I assume this feature 
of the results and define rup by 



1 _^ g+(fc+f)-g-(fc-f) 

777-7) rZ \ TP >^ TP 



(119) 



Lp 



^+F 



As mentioned above, vrip for the gapless cases merges 
smoothly with the rUp for the nea rby gapp ed states. The 
results are thus given by Eqns. (109-llC) for the range 
.75 < A < 2 for the ZDS potential. 



C. Why the free-field behavior? 

The fact that magnetic phenomena at T = can be 
described (to excellent accuracy) by free fermions of mass 
rUp needs to be properly understood and interpreted. 
For example one must resist the thought that perhaps 
by some further change of variables one could take the 
present hamiltonian and convert it to a free one. This is 
because if there were really an underlying free theory, it 
would be able to predict an activation mass vria and this 
would have to coincide with nip. We know within this 
theory, within Jain's approach, or from experiment, that 
these masses differ by at least a factor of two. 

I will now show that a single assumption about the 
form of the ground state energy, an assumption that is 
not equivalent to the free-field assumption or even to a 
quadratic dispersion relation in the gapless cases, will 
explain this behavior for gapped and gapless fractions. 
Consider E{S), the ground state energy as a function of 
S — nP. By rotational invariance it must have only even 
powers of S in its series. Assume the series is dominated 
by the first two terms: 



a 



E{S) = E{0) + 1^2 



(120) 



where a is the inverse linear static susceptibility. 

Consider first the gapless case. When dn particles go 
from spin-down to spin-up, 



dE = a S dS = a S {2 dn) 

h2 _ 1.2 

-(2 dn) 



Aw 



(121) 
(122) 



using the volumes of the Fermi seas. We see that dE 
has precisely the form of the kinetic energy difference of 
particles of mass mp given by 



1 



a 

TT 



(123) 



Thus rUp is essentially the static susceptibility, which 
happens to have dimensions of mass in d = 2. The state- 
ment that mp has no r-dependence in the gapped case or 
no spin dependence in the gapless case is the same as say- 
ing that the full nonlinear susceptibility does not depend 
on the spin S, which in turn means E(S) is quadratic in 
S. 

Note that the free-field form of dE comes from E ~ S^ 
and d = 2: in d = 3, we would have dE/dn ~ 5* ~ 
(k+F ~ '''-.f) which no one would interpret as a difference 
of kinetic energies. 

This general argument notwithstanding, it is worth 
explicitly considering the case in hand. 

First, one must not think that E{S) being quadratic 
in S as equivalent to assuming that the single-particle 
HF energies are quadratic in momenta. Consider the HF 
energies quoted earlier 



£ik 



±F) 



'/el) 



±F 



±F 



(124) 
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The quartic terms miraculously drop out in the energy 
cost of transferring a particle from the top of the spin- 
down sea to the top of the spin-up sea: 



2t* 77 

S = nP = n(l ) = -(p- 

P P 



2r) 



(132) 



It then follows that 



-a + . ^ - + b- 



(eVd) 



Kp 



kp 



(125) E{p -r,r)- E{p~r-l,r + l) 



K+F ''^-.F , A^+F k_p)[k^p + k_p) 

a ,^ 1- b — r-r^ (126) 



kp 



kp 






kp 



(127) 



using 



^+F 



k_p 



kp. 



(128) 



Note how d 



we would have k^p 



2 was essential to this argument: in d = 3 

l^_p — h,p. 

Thus the fc^ terms in £{k±) are not the cause of the S'^ 
term. However, a small k^ term in £{k±), corresponds to 
small quartic terms in E{S). 

To understand why the k^ term is so small, we turn 
to Eqn. ( |63| ) for Hq. Expanding the sin^ in a series, we 
find the fc° term is down by a factor of at least 15 (50) 
relative to the k^ term, at A = (A = 1), all the way up 
to k = kp. Presumably this feature (and its counterpart 
in the gapped case) persists in the HF approximation to 
H and keeps E{S) essentially quadratic, which in turn 
mimics free-field behavior. 

To really drive home the point, consider a problem 
where particles are free and have a dispersion relation 



S{k) = ak^ + bk^ 



(129) 



Let yn be density of spin-down particles and (1 — y)n 
particles that of spin-up particles. Thus 



S=[l- 2y)n. 
The total energy, as a function of y is 



(130) 



(>4:7Tny />4:7zn{l~y) 

E{y) ~ / dk^iak^ + bk'^) + / dk^{ak^ + fefc"*) 



a\y^ + {l-yf)+b\y' + {l-yf) 



(131) 



Note that the cubic terms in y cancel. Since y is lin- 
early related to 5', it follows E{S) is also quadratic in S. 
However a k^ term would have led to S"' terms in E{S). 

We wrap up this topic with one thought: even if CF are 
free or nearly so, there is no reason their kinetic energy 
should be quadratic in momentum. These particles owe 
their kinetic energy to electron-electron interactions, and 
given this fact, all we can say is that their energy must be 
an even function of k, starting out as fc^ at small k. What 
constitutes small k is an open question that is answered 
unambiguously here: our expression of the energy has 
substantial k'^ terms for momenta of interest. 

Let us now turn to the gapped case and verify that 
E{S) = E{Q) + f 52 implies that A(r) will be r- 
independent. First note that 



-{S^{p^r,r)~S\p-r-l,r + l)) 



{p-2r+p-{2r + 2)){p - 2r - p + 2r + 2) 
(2p - 4r + 2) 
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fn^l 


2 
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n^l 


a 


p2 



''-{p-2r-l)Apirf'f 



P 

We find that the r-dependence of Ap(r)'^'^^ drops out 
and gives 



A = 



2na 



If we write 

A„ 



eB* 



27m 



eB 



mp{2ps + 1) 
1 2?™ 



I'mp 2ps + 1 rripp 
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(134) 
(135) 



it implies, upon comparing to Eqn. (133) that, as in the 
gapless case 



1 



a 

TT 
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Analogously to the gapless case, we can show that if 
the HF energies vary with the LL index n as £{n) — 
an + bn"^, E{S) will be quadratic in 5*. Thus the CF-LL's 
do not have to be uniformly spaced for them to behave 
as if they were (with a spacing Ap) in T = magnetic 
transitions. 



D. Effective potentials for dirty systems 

Here I ask if it is possible that a ZDS potential with 
some effective A can describe the dirty system. First of 
all, I realize this cannot be true with respect to all observ- 
ables, if at all it is true for any. For example, if one were 
considering conductance, one knows the electron in a dis- 
ordered potential will typically get localized whereas no 
ZDS interaction will predict this. As for transport gaps, 
the present day samples, with a disorder broadening of 
the same order as the gaps, again preclude this possi- 
bility. Magnetic transitions, on the other hand, are con- 
trolled by total energies and one may expect that disorder 
will have a rather innocuous effect and can be represented 
in an average way by some translationally invariant in- 
teraction. I will show below that at least in some limiting 
case this is a reasonable approximation. I will show that 
even in cases where I cannot provide a similar argument, 
it seem to work. This needs to be fully understood. 
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As seen above, the critical fields are controlled by 
ground state energies in states of different polarization. 
Let us write the ground state energy as 

Eo = {n\H\n) (137) 

where \Q) is the ground state, and 

H=IJ -^,r (q) v{q)e-^^^^'l^ r (-q)- (138) 

Now add on a perturbation that couples the system to 
an external impurity potential $(g). To second order we 
find, upon disorder averaging. 



, or n or v. 



the 



E — E{) 



''\e-^ 



{2^Y 



., \{N\wmv^l{q) 

En — En 



(139) 



where N is any intermediate state and $g is the aver- 
age of $(g)$(— q) over realizations. We will now replace 
En — -Eq by a, the smallest excitation gap produced by 
ff. In our HF states, this is the gap to the next CF 
LL, ignoring the q dependence of the magnetoexciton, 
or the roton minimum if one includes it. This replace- 
ment overestimates the second order contribution. Since 
matrix elements to more distant LL are accompanied by 
higher powers of q, the error will be small if the disor- 
der potential (due to faraway impurities) has only long 
wavelength components. If we now use completeness, we 
see that the second order energy can be found by sand- 
wiching a hamiltonian, once again quadratic in ff, but 
with an effective potential 



iJeffiq) ^v{q) 



^Uq) 



(140) 



While the second term, due to disorder, need not be of the 
original form, what is important is that it is translation- 
ally invariant and makes a negative definite contribution 
to Veff ■ We could incorporate its effect by increasing A in 
the original ZDS term by a suitable amount. While such 
a replacement cannot reproduce all the effects of Ug// 
in detail, let us note that the energies depend on the 
potential v(q) principally via the corresponding Haldane 
pseudopotentialcZl Vi. (This is for 2s = 2. For 2s — 4, it 
is V3. In any event, one only dominates.) We may thus 
choose the effective A, so that it reproduces the dominant 
V corresponding to fe//- 

We see that at T = 0, in a state with a robust gap, to 
second order in the impurities, there is an effective A, if 
we ignore higher CF LL's. These are a lot of restrictions. 
Of these, only use of second order perturbation theory 
may by justifiable since the CF has a charge e* and not 
e, and correlations, which are taken into account from 
the outset, reduce the coupling to disorder by a factor 
l/(2ps -I- 1)^ at small q. The other restrictions cannot be 
justified at this point. I will however pursue the notion 
of an effective A for all fractions and at T > 0. 

Specifically, A will be extracted from one data point 
and used to explain the rest of the data from that sample. 
If the other data points differ only in the temperature T, 



the same A will be used. If it differs in , 
following scaling argument will be used.ca 

In a heterojunction, the donors of density n produce 
a confining linear potential of slope that goes as n. If 
one considers a variational wavefunction of the Frank- 
Howardo form ip{z) — A{w)zexp{—z/w) in the trans- 
verse direction, then the optimal w (to which A must 



be proportional), varies as w 
A = A/l varies as 



-1/3 



Consequently 



A 



-1/3^1/2 ^ 5l/6j,-l/3 



n'/^iy-'/\ 



(141) 



E. Comparison to data of Kukushkin et al . 

Kukushkin et al-H vary both n and B and drive the 
system through various transitions at T = (by extrap- 
olation) . The field B is always perpendicular to the sam- 
ple. We will compare the hamiltonian theory to these 
experiments by calculating the critical fields at which the 
ly — 1/2 and v — 1/4 systems saturate (P — 1) and the 
gapped fractions undergo transitions from one quantized 
value of P to the next. 

Let us recall that as far as these transitions go, the 
systems behave like free femions of mass rrip which is 
independent of the index r which labels how many LL's 
have reversed their spins in the gapped case or in the 
gapless cases, the size of the up and down Fermi circles. 
At and near the gapless states rrip may be fit by the 
expressions 



1 



m 



(2) 



m 



(4) 



= ^Cf(A) 






.087 
AV4 

.120 



(142) 
(143) 



where the superscripts on C refer to the number of vor- 
tices attached. 

I consider P'^'s at which the systems at 
1/4,2/5,3/7,4/9, and 1/2 lose full polarization (r == 
for gapped cases, saturation for the gapless cases) and, 
for 4/9, also the r = 1 transition, |3, 1) -^ |2, 2). 

I fit A to the v = 3/7 transition |3,0) -^ |2, 1) at 
B" = A.5T. 

Two points need to be mentioned in connection with 
the experiment. First, each of these transitions seems 
to take place via a narrow intermediate step with a po- 
larization half-way between the ones allowed by CF the- 
ory based-on spatially homogeneous states. Murthy has 
suggestedEy that these correspond to spatially inhomoge- 
neous states. I use the center of these narrow steps as the 
transition points for comparison to the present theory. 
Secondly, I use the actual values for Cp for these frac- 



tions, rather than the asymptotic value in Eqns. (142). 

I then obtain A3/7 ~ 1.42 on solving Eqn.(lll) which 
takes the following specific form here : 



B" 



2Ap[(3,0)-.(2,l)] 2(.0117) 



2me{e^/el) 



(eVd) 
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(144) 
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For transitions at other B± and n, I need the corre- 
sponding A's. One can argue, as per Eqn. (141) that 
A ~ I'^^/'^B^/^, from which it foUows that 



'^i^ — -^3/7 



B 

45 



1/6 



3 

7^ 



ll/3 ^1/6 



(145) 



Given A one finds B'^ using Eqn. (Ill) for gapped 
cases. 

For the gapless cases, there are two equivalent ap- 
proaches. First, at the critical field the Fermi energy 
of the up spins equals the Zeeman energy of the down 
spins: 
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(147) 
(148) 
(149) 



Using Eqn. (145) one solves for B'^ and obtains the values 
given in Table g. 

Equivalently we can write for the total ground state 
energy density E^{S), (where the superscript indicates 
that the Zeeman energy is included), 



E'{S) = ^S' 



B^S 



2m„ cos6' 



(150) 



where a — -k /nip. This expression is minimized (forP < 
1) to give P: 



P = 



cos 9 
19JB~X'^/^ 



cos 9 



!/=-, BinTesla (151) 

!/=-, BinTesla (152) 



Setting P = 1 gives the critical fields. 
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5.9 T 


2.62 


1/2 


saturation 


9.3 T 


11.8 T 


4.65 



interactions will increase the effective thickness with in- 
creased density. If I fit to the 2/5 point, I obtain similar 
numbers, with the agreement worsening as we move off 
in density from 2/5. Thus 3/7 was chosen as the fitting 
point since its density was somewhere in the middle of 
all the densities considered. 

An alternate approach is to attempt to calculate A 
ab initio using, for example, the Local Density Ap- 
proximatioUpO^DA) as done by of PMJ and numerous 
predecessorsE^. This method will however not include 
the effect of disorder. It typically gives a A that is half as 
big. We cannot attribute the entire difference to disorder. 
For one thing, the HF approximation tends to inflate A. 
Next LDA , as the name suggests, is an approximation. 
Lastly, LL mixing can account for some_Df the difference, 
though not much at these values of Ac3. It is possible, 
that due to all these, the effects of disorder are not that 
significant, for magnetic phenomena. The present strat- 
egy of determining A phenomenologically stands or falls 
depending on how well the fit to any one data point allows 
us to make predictions for other measurements made on 
the same sample. Tabic (^ suggests it is quite a useful 
point of view. 



VI. PHYSICS AT NONZERO TEMPERATURES 

r>o. 

So far we have seen the hamiltonian theory may be 
used to compute quantities such as gaps, particle-hole 
profiles, critical fields for magnetic transitions and so 
on. All such quantities have been readily computed using 
trial wavefunctions, giving numbers that are superior to 
ours. The main idea of this paper so far has been to get 
a 10-20% theory in which we see the underlying physics 
as transparently as possible and to resolve questions such 
as why CF behave like particles on some occasions. 

We turn to physics at finite T where the method 
has few rivals. Exact diagonalization is limited to very 
small systems and trial wavefunctions typically cover the 
ground state and very low excitations. The hamiltonian 
approach is able to yield, in the HF approximation, the 
polarization P and the relaxation rate 1/Ti for the gap- 
less states as a function of temperature. If A is treated 
as before (fit to one data point per sample) we will see 
it is possible to give a very satisfactory account of ex- 
periments up to about 1°K, which is of the order of the 
Fermi energy. 

The HF energy of a particle including the Zeeman en- 
ergy is the self-consistent solution to 



Table 6 Critical fields based on a fit at 3/7. The rows are 
ordered by the last column which measures density. 



Note that in rows above (below) 3/7, where I fit A, 
the predicted B'^'s are lower (higher) than the observed 
values, i.e., the actual A's are less (more) than what Eqn. 
(|145D gives. This is consistent with the expectation that 
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where the superscript on £^ reminds us it is the total 
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■gy including the Zeeman part, the Fermi functions 
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FIG. 8. Hartree Fock energies at i/ = 1/2 for up and down 
spins (upper and lower curves) at T = .3°K at B — 5.52 
and zero tilt. Note that they are not simply quadratic in 
momenta, and that at the chemical potential, indicated by 
the horizontal line, the two graphs have very different slopes, 
i.e., density of states. 



:(ifci) = 



cxp [(5|(fc) - ^)/fcT] + 1 



(153) 



depend on the energies £±{k) and the chemical potential 
/i. At each T, one must choose a fj,, solve for £±{k) till a 
self-consistent answer with the right total particle density 
n is obtained. From this one may obtain the polarization 
by taking the difference of up and down densities. As 
usual we use the ZDS potential for which 
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The computation of 1/Ti is more involved. The ques- 
tion we ask is the following. The fermions are in a quan- 
tum well, with their density varying across the width. 
So the nuclear relaxation rate will be a function of po- 
sition. Consider a nucleus at the center of the quantum 
well, (as well as the x — yplane) where the density is the 
largest. Let us call this point the origin and let l/Ti be 
the relaxation rate here. The theory predicts 



1 / JV'™"^ 

— = ^TxksT -^ — 
Ti V ^ 

dE {'^-^^4S^ 1 p+{E)p^{E)F{k+,k^) (155) 
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where Eq is the lowest possible energy for up spin 
fermions, and iT™"^ is the measured maximum Knight 
shift for the fraction i/ = 1/2 or 1/4. 

Here is a rough description of the derivation, the de- 
tails of which may be found in the Appendix 6. Suppose 
for a moment we were dealing with electrons and not 
CF's. The Knight shift at the chosen point, the origin, 
will be determined by the spin density there. The same 
parameter enters the l/Ti calculation quadratically. This 
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FIG. 9. Comparison to the work of Dementyev el al. The 
value of A is fit to P at 300 mK, B± = 5.52 T and the rest 
follows from the theory. Notice the correlation between the 
curvature of l/Ti and the limit of P as T ^ 0°^". 



is why X™"^ enters the answer. The idea is that KJJ^"-^ 
is not calculated ab initio but taken from the same ex- 
periment. The density of states and Fermi factor are 
standard. The only new feature here is the presence of 
F{k+,k-) which reflects the fact that the spin density 
has to be projected into the LLL when going to the CF 
basis. The effect of this factor (which is none other than 
the e~' ' which appeared on the projected charge den- 
sity) is to suppress processes with momenta much larger 
than l/l, as these have no place within the LLL. 



A. Comparison to experiment 

We now compare to some experimentSpat i^ — 1/2 and 
T > 0. Consider first Dementyev et at3. From their 
data point P = .75 for B = Bj_ ^ 5.52T at 300 mK I 
deduce 



A== 1.75. 



(158) 



I have once again chosen instead to match my HF results 
with the above data point, (which gave A = 1.75) and 
see to what extent a sole parameter A, can describe P 
and 1/Ti for the given sample at a given B±, but various 
temperatures and tilts. 

Since there does not exist a model, including disorder, 
that describes how A should vary with tilt I include no 
such variation. 

Dementyev et al find K^^l'' = 4.856 • lO"'^"^, which is 
believed to describe a saturated system at P = 1. They 
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FIG. 10. Comparison to Melinte et al. with A fit to P at 
60 vnK, B± = 7.1 T. Some typical error bars are shown. 



estimate that K^,"^^ 



.953i^]^;^^, 



which is what we need 
here. It is assume d tha t the nuclear wavefunction u_b(0) 
(defined in Eqn. (224) Appendix 6) does not vary be- 
tween v = 1/3 and v = 1/2. Given this information, 
1/Ti follows. 

The top and bottom halves of Figure compare the 
HF calculation of 1/Ti and P respectively, to the data. 
(The graphs for 1/Ti differ slightly from those in Ref.B 
since the present calculation treats the spin of the CF 
more carefully. The 1/Ti graph at 5.52 T appears a little 
jagged since it was computed at just six points which 
were then connected. This is not apparent in the tilted 
case since the points lie on a straight line.) 

Dementyev et al had pointed out that a two parameter 
fit (using a mass m and interaction J), led to disjoint 
pairs of values for these curves. Given that H is neither 
free nor of the standard form {p^ /2m + V{x)) this is to 
be expected. By contrast, a single A is able to describe 
the data here rather well since H has the right functional 
form. Given how the theory fits the polarization data up 
to the Fermi energy of ~ 1°K, it is clear that changing 
the data point used to fix A will be inconsequential. 

If P were computed from the LDA value A ~ 1, it 
would be down by 15-50 % as T drops from 1°K to 0°K. 
The present work establishes a phenomenological, non- 
trivial and nonobvious fact that a single A parameter, 
(like g or e) determined from one data point, can de- 
scribe both P and 1/Ti for the given sample under a 
variety of conditions. That the fitted A is larger than the 
LDA value makes sense, as both disorder and LL mixing 
will lower the gap and raise A. As mentioned earlier, it is 
not clear exactly how strong the disorder corrections are 
since there is LL mixing and some intrinsic errors in the 
LDA. 

Consider next sample M280 of Melinte et aO which 
had P = .76 at .06°K and B ^ B± ^ 7.1T, from 
which I deduced A = 1.6. Figure ( [lO| ) compares my T 
-dependence with data. The initial rise of P with temp|ei;- 
ature was also seen by Chakroborthy and PietlianeanEJ. 

Note that agreement is quite poor for the tilted case. 
There are confusing aspects of both the theory an exper- 



iment. In the experiment, one may ask why the polar- 
ization does not increase with increased tilt and hence 
increased Zeeman coupling. Of course this will happen 
if P = 1 to begin with, but it is not, it is clearly below 
unity, say 80%. On the other hand, the theory for tilted 
fields in not in good shape either. First of all, orbital 
effects have to be considered due to the tilt. The thick- 
ness parameter A can be affected by it. As pointed out 
by Jungwirth in a private communication, once there is 
an in-plane component of P, the problem is no longer 
rotationally invariant. This means that our states are no 
longer HF states and get scattered into each other by the 
potential. Presumably these effects were negligible in the 
case of Dementyev et al but not at nearly twice the tilt 
of 61°. No attempt is made here to take into account all 
the effects of the tilt. Instead I include just the increased 
Zeeman coupling and hope for the best. 

For the benefit of others who measure l/Ti at i/ = 
1/2 in the future on similar samples, I give some very 
approximate formulae (to be used for zero or small tilts). 
From Figure (||), we note that in general, the graphs of 
l/Ti become linear and parallel for temperatures above 
.3 °K. In this region we can write 



dT 



■ IQ-^s-^ \°K] 



with K the Knight shift in KHz and 



lOio/cm^ 



for T > XK 
(159) 

(160) 



(In this appro ximate f ormula, I ignore the A depen- 
dence of Eqns. (155-157), and the distinction between 
the average and maximum Knight shift.) 

The graphs do not generally obey the Korringa-like 
law because as T — > they are sublinear (super linear) for 
saturated (unsaturated) cases. Only the critical case with 
P(0) ^ 1 as T ^ is hnear. For T > .3 ° iC, (which 
in general must be replaced by either the energy gap or 
energy overlap between the up and down Fermi energies) 
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K 

n 



- = [3 Ti^K) + C] 
J-i 

C = (critical) 

= > (unsaturated) 
= < (saturated) 



■ 10-^s^i 



(161) 

(162) 
(163) 
(164) 



For the critical case (only) , we have Korringa law 

>T 2 



1 



Ti T 



= [3 T{°K)] 



lO-^s-^{° K)-^ 



(165) 



For Dementyev et al C ~ 1. This value may be used 
as a first approximation, for example, if Melinte et al 
measure l/Ti on sample M280. F or more accurate results 
they must solve Eqns. (155-157) with A = 1.6. 

The present formalism has been applied by G. Murthy 
to calculate the_T-dependence of polarization in the 1/3 
and 2/5 statesH 
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VII. CONCLUSION 

This paper described in detail a formalism for describ- 
ing FQHE states, in which the hamiltonian and various 
operators were expressed in terms of the quasiparticles, 
the Composite Fermions. It was shown that the formal- 
ism could be used to calculate a variety of quantities at 
zero and nonzero temperatures to some reasonable accu- 
racy, and to resolve matters of principle (some of which 
could not even be posed otherwise), such as the internal 
structure of the CF, the reason it sometimes appears to 
be free when it is not, and how theory is to be compared 
to experiment in the quantum Hall problem. 

First a review of the wavefunction approach to CF was 
given. It was seen that whereas in the Laughlin fractions, 
CF was evidently an electron bound to 2s vortices, the 
situation was more complicated for the Jain series due to 
the action of the LLL projector V. This caused the vor- 
tices to move off the electrons and also got rid of many 
of them. It was not clear how the vortices paired off with 
the electrons, say at i^ = 2/5, when there were 1.5 (non- 
Pauli) vortices per electron in the projected state. Since 
the charge of the CF was the same (e* = l/(2ps + l)) be- 
fore and after projection, whatever paired with the elec- 
tron still had the charge of 2s vortices, though it could 
not be associated with zeros of the wavefunction in any 
simple way. What is this object and how is one to incor- 
porate it in the theory? 

The hamiltonian theory was seen to provides the an- 
swer. In this theory we enlarged the Hilbert space to in- 
clude additional degrees of freedom, accompanied by con- 
straints. These new degrees of freedom (prevented from 
having any density fluctuations by constraints) turned 
out to have charge e^ = —2ps/{2ps + 1) and paired with 
electrons. They were referred to as vortices, for want of 
a better name, but in view of what was said above, they 
are not related in any simple way to zeros of ^ lll, the 
electronic wavefunction in the LLL. 

To solve the theory with constraints, the HF approx- 
imation was used. Free particle and hole states of CF's 
in a reduced field B* = B/{2ps + 1) were seen to be HF 
states. All matrix elements were evaluated in closed form 
and usually a single numerical integral gave the numbers. 
As for the constraints, the method proposed with Murthy 
was used: the charge density operator was written as a 
judicious combination of the transformed and LLL pro- 
jected electronic charge density and the constraint, which 
is just the vortex charge density, the judicious combi- 
nation being the unique one (as q ^ 0) that obeyed 
Kohn's theorem. When this combination was expanded 
in q, the monopole term was e* and the next term cor- 
responded to a dipole moment Pz x k. The operator 
approach thus gave a precise meaning to the CF dipole 
moment, in terms of how it coupled to an external poten- 
tial. This coupling did not change when second quantized 
with Fermi fields. Since this procedure took into account 
the most important effects of the constraints (away from 
T = and ultralow frequencies w ~ q^), constraints were 
neglected thereafter. 

In Section IV activation gaps for fully polarized states 



were calculated for the Zhang Das Sarma (ZDS) potential 
as a function of the parameter A. Analytic expression 
were derived for all fractions of the form p/{2ps + 1) and 
compared to the Monte-Carlo work of PM J based on trial 
wavefunctions for s — l,p — 1,2,3,4. It was found that 
the numbers were within 10-20% for A > 1. j— . 

A comparison was made to experiments of Du et a£3. 
Rather than try to compute A ab initio, it was fitted to 
the three gaps at 1/3, 2/5, /3/7 and was seen to take 
the values 2.4, 2.4, and 2.0. 

It was found that normalized activation masses (scaled 
by 1/vB) were not too different for s = l,s = 2, i.e., 
near v — 1/2 and v = l/4r,in reasonable accord with 
the experiments of Pan et aO.It was pointed out that no 
deep reasons existed for their exact equality. 

For the gaussian potential the numbers were in excel- 
lent agreement with PMJ except for v = 1/3. In was 
clear that although the extended theory was defined for 
all length scales, it gave good numbers only for soft po- 
tentials, ones that were smooth within a magnetic length. 

A comparison to the exact diagonalization work of 
Morf et al was made for fractions 2/5 and 3/7. The HF 
results of this theory lay consistently above their num- 
bers. The differences were somewhat larger than in the 
case of PMJ. 

In the hamiltonian formalism we can make precise the 
question of whether or not CF are interacting since H 
naturally separates into a free and interacting parts Hq 
and Hi. It was seen that turning off Hj made a sizeable 
difference to activation gaps. 

The theory was used to compute the profiles of charge 
densities in states with a single particle or hole and com- 
pared to the unpublished work of Park and Jain. It was 
seen that while the salient features were reproduced with 
ups and downs at the right places, the amplitudes were 
not as pronounced. The best fraction was 2/5, and not 
surprisingly, this was also the one where the gaps came 
out best. The case i/ = 1/2 and the dipole moment of 
the CF were analyzed in some depth. 

It was found that the activation gaps for the fully po- 
larized states obeyed particle-hole symmetry to an excel- 
lent approximation. This was a nontrivial result since 
the corresponding CF states were very different as were 
the expressions for gaps. 

Polarization phenomena at T = were the subject 
of Section V. By varying the density and field at fixed 
filling, or by placing the sample in a tilted field it is pos- 
sible to increase the Zeeman coupling and drive the sys- 
tem though many magnetic transitions. For the fraction 
p/{2ps + 1), CF theory has states with p — r LL of spin 
up and r levels with spin down, r = being the fully 
polarized case. Critical fields B'^ at which the system 
would jump from one value of r to the next were com- 
puted. It was found, as Park and Jain did, that it was 
possible to fit all the numbers very well by assuming that 
CF were free and occupied LL with a polarization gap Ap 
or a corresponding polarization mass mp. Yet we know 
CF cannot be free, given that the activation mass and 
gap are substantially different from these values. It was 
shown that rotational invariance and d = 2 conspired to 



24 



create this impression of free fermions. In particular, for 
the gapless case, while the individual HF energies of spin 
up and down fermions had sizeable fc^ terms in the for- 
mula for £±{k), the energy cost of transferring a fermion 
from the top of one sea to the top of the other equaled 
that of free fermions with a quadratic dispersion k'^ /2mp. 

The theory h«is compared to the experiments of 
Kukushkin et aJ13 who varied both n and B± to drive 
the magnetic transitions. Some plausibility arguments 
were given for why an effective, translationally invariant 
potential might be able to describe magnetic phenomena 
under restricted conditions. It however used to describe 
all magnetic phenomena: an effective A was deduced from 
one data point (here one transition), and used to predict 
all the others, using scaling laws. The numbers agreed 
to within 10 — 20%. While the arguments for an effective 
potential description were trustworthy only in a limited 
region, the procedure seemed to work in a wider region, 
including gapless states and T > 0. This deserves to be 
understood. 

The most important results of this theory, not limited 
by finite size, involve the computation of relaxation rates 
and polarizations as a function of temperature. These 
were discussed in Section VI. Finite T HF equations were 
derive analytically and solved numerically to yield these 
quantities. They were compared to the work of Demen- 
tyev et al and Melinte et al. In the first case, the mea- 
sured polarization at 300 mK for a perpendicular field of 
5.52T was used to fix A = 1.75. Using this value the po- 
larization and relaxation rate 1/Ti were computed for a 
range of temperatures going up to 2°K, and 1°K respec- 
tively, where the latter is roughly the Fermi temperature. 
The agreement was very satisfactory. For a tilt of 38.3°, 
the agreement was again good till about 1°K. In their 
paper Dementyev at al had pointed out that fitting these 
four data sets to a model with a mass term m and Stoner 
enhancement J lead to four disjoint islands in the param- 
eter space. In the present case a single A = 1.75 seemed 
to describe all four sets quite well. It was claimed this 
was because the H used was of the right functional form, 
with an unusual kinetic and interaction terms chock full 
of momenta and currents. For any one graph it may be 
mimicked by a hamiltonian of the (to, J) form, but these 
numbers would then vary from set to set. In particular, 
to get the right relaxation rates, with just one param- 
eter, A, one needs the unusual dispersion relations (far 
from quadratic in momenta) and the corresponding un- 
usual density of states that arise naturally herej— . 

As for the polarization data of Melinte et all3, there 
was reasonable agreement for the untilted sample, but 
not the tilted (by over 60°) case. The latter could be 
attributed to theory, which is not designed to handle such 
large tilts, or to the data which exhibit some unusual 
features discussed in the text, or both. 

An approximate Korringa-like law for states that are 
just fully polarized at T = and the slope of l/Ti versus 
T for the general case for T > .3 °K were provided to 
allow future experimenters (working on samples similar 
to that used by Dementyev et al) to make a quick and 
approximate comparison to this theory without having 



to solve the finite-T HF equations. 

In summary it was shown that the hamiltonian theory 
of CF provides a comprehensive and analytical scheme for 
describing the low energy physics of the FQHE states, 
clarifying concepts regarding the internal structure of 
CF's and how they are to be coupled to external poten- 
tials, and computing a variety of quantities at zero and 
nonzero temperatures to an accuracy of 10 — 20% un- 
der typical conditions, and sometimes considerably bet- 
ter. The deviations could be due to the HF approxima- 
tion or to the hamiltonian itself. In comparing to mag- 
netic experiments, if one extracted a single parameter A 
that characterized electron-electron interaction in a given 
sample from one data point, the theory gave a reasonable 
account of other data from that sample. 



A. General philosophy 

I conclude with some remarks that put the present 
approach in perspective and underscore the differences 
between the FQHE problem and others. Let us begin 
with the fact that in the FQHE restricted to the LLL, 
and in the absence of disorder, every physical quantity 
is a functional of the electrostatic interaction between 
electrons, since the kinetic energy is quenched. I will 
address the important question of disorder shortly. If this 
interaction is pure coulomb, it is a zero-parameter theory. 
If it is modeled by the ZDS potential as done here, all 
observables - nia, nip, 5°, P{T), 1/Ti - are functions of A. 
Once a single data point is known, A may be determined, 
and from it, all others calculated (using scaling laws if 
needed). While I chose the ZDS potential to make this 
poinjt. one could use another, say the one used by Morf 
et alB. The the numbers in question are predominantly 
controlled by one of Haldane's V's, and these are different 
ways of varying it. 

Why do we not do this all the time? Consider for ex- 
ample a Fermi liquid. Why do we bother with Landau's 
F functions? Why do we not start with the coulomb in- 
teraction and calculate everything with zero parameters? 
Or if we wanted to model it with a short range force, why 
don't we employ the Hubbard model, calculate any one 
observable, fit it to data, extract the Hubbard U and 
predict everything else? The reason is well known: we 
do not know how to go from the model to the physical 
quantities except in perturbation theory. Let us restate 
this in Landau's language. Suppose we turn off U. The 
states, at least the low lying ones, are labeled by free field 
theory. We now turn on the interaction and the states 
and energies evolve to those of the interacting theory. 
We finally should end up with the quasiparticle basis, at 
least near the Fermi surface. Unfortunately this change 
of basis from free to interacting theory cannot be carried 
out in practice except perturbatively. 

The same thing happens in S-matrix theory, which 
appears to have guided Landau's intuition. The in and 
out states carry free particle labels and their dot product 
gives the elements of the 5'-matrix. They are obtained by 
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taking particles that are infinitely separated (and nonin- 
teracting) in the distant past or future and evolving them 
(forwards or backwards) in time to t ~ 0. Unfortunately 
this evolution too can be done only pcrturbatively. 

In the FQHE problem the situation is a lot better, 
though it looks bad to start with. The usual idea of 
starting with free electrons and turning on interactions 
is doomed from the start since for v < 1, the free elec- 
trons do not have a unique ground state. But lurking 
in the background are again some free states or wave- 
functions. These are the CF states with p filled LL and 
some low energy excitations. They get mapped to the 
electronic wavefunctions in Jain's approach by a trans- 
formation that consists of attaching 2s-fold vortices with 
the Jastrow factor J {2s) and projecting with V. The 
resulting wavefunctions arc known to be excellent, with 
nearly unit overlap with the results of exact diagonal- 
ization. If we know the electron-electron potential (in 
the absence of disorder) we can calculate anything re- 
lated to the ground state and low energy excitations us- 
ing these. In the hamiltonian approach we go the other 
way: from electrons to CF. We first do 2s-fold flux at- 
tachment by the CS transformations and this becomes 
vortex attachment following the additional transforma- 
tions Murthy and I developed. The states not only carry 
free labels, they are simply free (in the HF calculations); 
it is the operators that become complicated. They are 
however simple enough to work with analytically, and 
lead to less accurate numbers than Jains 's approach. 

All this is possible because of the one major difference 
compared to the Fermi liquid and the 5'-matrix prob- 
lems: the transformation here is discrete and given by 
the attachment of 2s vortices (before projection, which 
leads to a complicated wavefunction, but is fully imple- 
mentable). The CF know only about p, the transforma- 
tion about 2s, and the electrons know about p and 2s. 
It is this discreteness, (absent in the Fermi liquid or 5*- 
matrix problem), that eventually works in our favor, that 
sufficiently constrain the problem and lead us to the right 
answer. For example, in Jain's approach analyticity and 
Fermi statistics fix the vortex number to be 2s and in the 
hamiltonian approach the same is true as well, and addi- 
tionally, the LLL algebra pointed to the exponentiation 
which was the nonperturbative step that had to exist in 
any transformation linking electron to CF. 

There is a mistaken belief that CF have to be free to 
be useful. Note that Jain is able to calculate a slew of 
quantities with his wavefunctions and to a lesser accu- 
racy, so can the hamiltonian method. Neither method 
requires CF to be free. What is free about the CF is the 
label for the states. While this is true in S-matrix theory 
and in Landau theory, what is special here, as explained 
above, is that the transformation from the free states to 
the interacting counterparts is known to a great accuracy. 

I believe this happy situation, of being able to go from 
the free to the interacting states in the FQHE, deserves 
further exploitation in comparing to experiment. If there 
were no disorder, what we should try to extract from the 
data is the electron-electron interaction, v{q), which is a 
c-number function that multiplies the operator quadratic 



in p^ and defines H . By extracting v{q) I mean some pa- 
rameter like A of the ZDS potential or b in the Morf et al 
case. As mentioned earlier, the physics is dominated by 
one Haldane pseudopotcntial, {Vi for 2s = 2) and these 
parameters just control it. It docs not seem fruitful to 
fit the complicated CF hamiltonian to a string of stan- 
dard or nonstandard operators. The hamiltonian for the 
CF's is not of the canonical form with a quadratic kinetic 
energy and some interaction terms of the usual (density- 
density) form. There is no reason such should be the 
case, given the internal structure of CF's. Rather, the 
free and interacting parts of the CF hamiltonian form a 
monolith, of unusual functional form, fully determined 
by the density operator written in the CF basis. While 
it is possible to fit any one experiment with one set of 
standard interaction parameters, (like m and J that De- 
mentyev et al used) the fact that the functional form is 
wrong will manifest itself, as they found out, in the need 
for many disjoint sets of parameter for different measure- 
ments. (This is analogous to the fact that if we try to fit 
the kinetic energy of a relativistic particle as a function 
of its velocity to a nonrelativistic form, we will need a 
velocity dependent mass, while if we fit it to the correct 
form we will extract a fixed rest mass.) Likewise, if one 
extracts A (or its counterpart for another potential) from 
the data, a single simple result is more likely to emerge. 
This A will of course have the usual variance of 10 — 20% 
characteristic of this theory. 

Now it is time to wake up and .smell the disorder. 
Does it completely destroy the approach presented here? 
Is it possible that even the disordered system can be mod- 
eled by a pure potential with some effective A or 6? The 
answer is clearly negative if one wants to describe every- 
thing. For example no translationally invariant interac- 
tion will predict localization. As for activation gaps, the 
data suggest that disorder effects are quite strong: the 
disorder broadening (measured by the negative intercept 
of the gaps versus effective field B*) is comparable to the 
transport gaps. As B* -^ 0, the pure system gaps never 
vanish while in experiments they do, implying no effective 
potential exists. If at some future date, we get samples 
with even lower disorder, we may be able describe these 
gaps with an effective potential, even for small i?e// • 

For the present, for magnetic transitions at T = 0, 
which depend only on total energies, I have tried to ask if 
an effective potential exists which can subsume the effects 
of disorder. I gave crude arguments that showed that in 
the limit of weak and smooth disorder one could get the 
disorder averaged ground state energy to second order by 
replacing w(g) with a v^ff (g). While Veff did not have the 
same functional form as v, it was always weaker, which 
meant the effective A (which produced the corresponding 
dominant pseudopotential V) is larger. / have explored 
the possibility, of determining an effective A from one 
data point and using it (combined with a scaling law if 
necessary) to predict the other data on that sample, for all 
magnetic phenomena. While the plausibility arguments 
were given for T = and that too under a variety of 
restrictions, I have tried it for all magnetic phenomena. 
What I find that if one is prepared to work with a 10 — 
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20% theory, such a program actuaUy works. If such an 
effective potential is found from static or equal time data, 
it can be imported to the hamiltonian scheme to do other 
things such as to dependent response functions and T > 
physics. 
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APPENDIX 



1. Matrix elements 

Many of the calculations performed in this paper deal 
with the preferred density p^. In second quantization we 
write it as 

P \W = / , '^m2n2 ™ini/'m2n2:niini (,ioDJ 

where d\^ „ creates a particle in the state \m2 112) where 
m is the angular momentum and n is the LL index. They 
are related to the CF cyclotron and guiding center coor- 
dinates, R and r/ as follows. Let 
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where I* = //VT^c? is the CF magnetic length. These 
obey the oscillator algebra 
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Similarly we define, in terms of the cyclotron coordinates. 
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which obey the oscillator algebra 

[a,at]=l (171) 

given 

[r^x,rjy]=tl*^. (172) 

The states \Tnn) are just the tensor products 
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where |00) is annihilated by both a and b. 

I will now evaluate matrix elements of e^'^^'^ and 
g-jq-77 g^^^ show that 
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where 
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L is the associated Laguerre polynomial, and mi > ni2- 
If ?7ii < 7712 one may invoke the relation 
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To establish Eqn. (174), consider the coherent states 
|.).e^^^|0) = ^M.,"^ (177) 
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with the inner product 
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First we write from the definitions given above 



(^le-^'i-^U) 



00 CXD 

E E- 

7711— ^2—0 

i?(z,z,q). 



=m2 ^mi 



(7772|e-''l-^|777i) (179) 

(180) 



On the other hand 
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Comparing Eqns. ( 179[ - 183| ) and matching powers of z°- 
we obtain Eqn. (174) if we recall 
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To establish 
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where |p) stands for the (ground) state with p-fiUcd LL, 
and i, / label single-particle excitations on top of this 
ground state. We want to show that this matrix element 
vanishes if i 7^ /, i.e., the hamiltonian does not mix these 
putative HF particle states. (This result was established 
for the small q theory by G. MurthyJla.) The proof, which 
relies on just the rotational invariance of the potential, 
applies with trivial modifications to the hole states, i.e., 
to 



(again for tt-i > 77,2) we just need to remember that the 
commutation rules of the components of ry have a minus 
sign relative to those of R, which exchanges the roles of 
creation and destruction operators and hence g+ and g_ . 
Now we consider matrix elements of p, x, ff . As 
a first step, let us express the operators Re and Ri, in 
terms of CF guiding center and vortex coordinates R and 
rf. We have seen that in the CF representation 
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This matrix clement in question takes the schematic 
form 
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if we recall P — P^(l — (?). 
It can similarly be shown that 

,0 z X n* „ , 

Thus in first quantization 
P^Yl exp(-7q • Ri) exp(-iq • rj^c) 

i 

X = X! exp(-iq • RJ cxp(-iq-?7,/c) 



(187) 



where 1 stands for r77i?7,i and so on, and J stands for an 
integral over a rotationally invariant measure: 



1 f d'^q -a^l^/2 



(196) 



(188) 



(189) 
(190) 

(191) 



/ = exp 



/ -q'l 



2;*2 



\^ps(2ps + l) 



(Vortex form factor) (192) 



Now we use Wick's theorem and perform pairwise con- 
tractions on the vacuum expectation value, bearing mind 
that 

• We cannot contract the indices 1 and 2 or 3 and 4 
since this will require that g = at which point the 
measure (which contains the potential) vanishes. 

• If we contract i and / we already have the desired 
result. 

Here is a representative of the contractions we can get: 
Pi2P34<5/i(l - n()523{l - n^)K{l ~ nf ). (197) 



Armed with Eqns. (Il74| and |l86| ) we may finally write ^^^^^ ^f ^^ ^^^ pg^^j function for the LL labeled by r7i 



-^7712^2 ;mini 



77f =0(p-l-77l) 



(198) 



]j}2^„^x/2 (zis+L 



7771 ! 



H, \ mi-m2 



V2 



Tmi-m2( 
m2 



772! f —icq^r 



* \ ni— TI2 



"1! V V2 



-c'-f- 



—iq^r 



V2c 

rm2'nii V/V rn2ni 



ix) 



and so on. Since f — I and 7 — 4, the (1 — 77^)(1 — 77!") 
1. The integrand assumes the form 



00 00 

E E 

m2=0 n2=p 



P/2(q)p2j(-q) 



E /'"/»2(q)Pn2n,(-q) 
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Superscripts on p™^^^ and p^^^^ which will be apparent 
from the subscripts, will usually be suppressed. 



2. Proof of Hartree-Fock nature of trial states 



Consider 



/ , Pmjm2V^)Pm2mi\ ^.) 

.7712—0 



where we have also used the fact that e~^'^'^ ■ e~*'i'^ = / 
in doing the sum over 7772, and -Fd?]) is some rotation- 
ally invariant function. It follows that every term in the 
sum over 772 vanishes unless 77/ — Ui due to the angular 
integral in q. 
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3. Particle-hole profiles 

Here we deal with both gapped and gapless states. 



Let us create a particle in the p+ 1-th LL and ask what 
the corresponding charge density looks like. In q space 
this is given by 



(r(q)) = E(pi^^^l^24iP)pi2(q) 



(199) 



1.2 



where /i labels both n = p and to = of the created 
particle, and 1 and 2 stand for the double labels summed 
over in the definition of ff. 

Now we do the contractions bearing in mind that we 
should not contract 1 with 2 since this gives the back- 
ground charge. This gives 



{r{ci)) = J2('^A){d2di)pi2 



1,2 
1,2 

= Pfi^i 

,-,2;*2„2 
'qH*^ 






-c^fe-"' /*- Lp 



2c2 



(200) 

(201) 

(202) 
(203) 

(204) 
(205) 



It is straightforward to go to real space by Fourier trans- 
form. The figures depict the charge so obtained plus the 
background charge v / {2i:P) in units of the background 
charge. The result is 



where 

ny) = 



1 



1 



vdve-y"/^Uvr/l)F{y) (206) 



g-y^(2ps+l)/4 / -t,2c"(2ps+l)/4^ 



y^c^(2ps + l) 



_^2f^-y\2ps+l)/4c-^^ 



r-,,2 



r(2ps + i) 



2c2 



The computation of the hole charge is analogous. To 
facilitate comparison to Park and Jain I have placed the 
particle and hole at antipodal points on a sphere of the 
same radius as the one they used. 



b. The CF atv^ 1/2 

As in the gapped case, we now place a fermion just 
above the Fermi sea. As explained in the main body 



of the text, we need a superposition of states of fixed 
Py — Po a-iid a sum over all values of Px to localize the 
CF at a; = 0. Since the dipole always has a nonzero x- 
component in this superposition, it does not get washed 
out. 

For this case, with A* = 0, the LL expression for 
charge simplify greatly : 



^''(q) 



4^ 



(-2^)sin(l^)4_^dk 



(207) 



In other words, the pair of labels (m,n) is replaced by 
a momentum vector q. The density operator (for each 
particle) 



pP(q) = g-ii<i-'^+(lxP ' /2) _ g-''(q-r-qxp 1^/2) 



(208) 



only connects states differing by q. Carrying out the 
Wick contractions we get to the following result (drop- 
ping irrelevant constants) 

{pix,y)) ~ /"dp,dfee^«-"e-«'''/4 (e-^^^^''''/^ _ ^^9.^0/^2^ 



exp 



x-poP/2 



-\ 2' 



P 



- [Po -^ -Po] 



(209) 



Similar methods may be employed to calculate the 
structure factor S{q). 



4. Activation gaps 

Now we need to find the energy cost of producing a 
widely separated particle-hole (PH) pair. This will be 
done by evaluating 

A„ = (p + P\H\p + P) + (p + H\H\p + H) 

-2(p|i?|p) (210) 

(211) 



= j E(P) + E{H). 

Jq 



where P denotes a particle added to the state labeled 
/i = (n = p, jTi = 0) and H denotes a state in which a 
hole has been made in the state p = {n = p — l,m = Q). 
Let us consider 



E{P) = {p\dA^d2dld4dl\p)pi2P34 
In performing the contractions we 



(212) 



• Do not make any contractions within H. This gets 
rid of £'0 ~ (p|_ff |p), the ground state energy. 

• Do not contract 1 with 2 or 3 with 4 since v(Q) = 0. 
We end up with 

Sf,3Si4S2i,{l - nf )(nf')(l - nf )] pupsi 
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Since A — ii — 1 in the first term, we can drop (1 — 
nf )(1 — nf) and for similar reasons (1 — 77,3 )(1 — rij ) in 
the second giving us 



00 00 



^(^) = E E PM2(q)p2M(-q) 

7712 — 712— P 

00 p—l 

- E E P2M(q)PA'2(-q). 

7772—0 TI2— 



(213) 



(214) 



Since the sum over m2 is unrestricted, we can use com- 
pleteness and e~*i'^-e"*'^'^ = / to get rid of the m-index 
altogether. Thus we end up with 



p-i 



E{P)^(J2\Ppn\'-Y.\PP- 



\n— p 



n=0 



p-1 



{p\p{q)pi-q)\p)-2j2\ppn\' 



n=0 



(215) 



(216) 



A similar calculation for the hole state gives (upon 
dropping the ground state energy as usual) 



EiH) 



where 



p-i 



-{p - l\piq)p{-q)\p - 1) + 2 ^ Ipp^i, 



n=0 



i\p{q)pi-q)W) = E l^(*)'^ 



(217) 



(218) 



Putting all the pieces together, and recalling the various 
matrix elements, we obtain Eqns.BftJTq). 



5. Critical fields for magnetic transitions 

We need to calculate 

Eip - r, r) = (p - r, r\H\ p - r, r) (219) 

the energy in a state with p — r spin-up LL's and r spin- 
down LL's. Since the HF calculation for the spinless case 
is very similar, this treatment will be brief. We write 



H =^ d\d2dldiPi2P34: 

1234-^9 



(220) 



with the understanding that a label like 1 stands for the 
triplet (nijTOi, Si), s being the spin. The matrix ele- 
ments pij are defined by 



P12 



(He- 



and as a result 



-iq-R. 



jq.jyc _ ^2^g-*q-?7/cl |2^ 



' Pn\r. 



E{p — r^r) = 

I ^ nf(.)(l-n2^(.))|p„,„j2^(m|/|m) 

^ nin2S m 

'• V ' 

=n/p 

where we acknowledge the fact that the occupation fac- 
tors nf and n|" can depend on the spin. We have also 
used the fact that the sum over all values of m is the 
degeneracy of each CF-LL, n/p. Carrying out the sums 
over ni and n2, we obtain 



E(jp — r,r) — 

p—r — l 

E i^'L\piq)p{-q)\ni) 



n 

Pjq 
r-1 



ni=0 



p—r — 1 

/ , \Pnin2 

"1,"2=0 



r-1 



^{ni\p{q)p{-q)\ni) - ^ \pn,n2 

ni,n2— 



ni=0 



It is now straightforward to compute the critical field 
for the transition |p — r, r) -^ |p — r — 1, r -I- 1) by in- 
voking 



E{p — r,r) — E{p — r — l,r+l)=g 



2?Tip 



B''- (221) 
P 



6. Calculation of 1/Ti 

As explained in the main text, our strategy for com- 
puting 1/Ti will be to compute K™""^, the Knight shift at 
the center of the well, in a fully polarized sample, in terms 
of an unknown nuclear matrix element squared |u(0)p, 
and then to express l/Ti (which depends on |u(0)| ) in 
terms of the measured value of i^™"^ . 



a. Calculation of Knight shift 

In first quantization, the hyperfine interaction of elec- 
trons with a nucleus at the origin is 



Hhf^Y^nleS-IS^Q) 



(222) 



where 7e and 7„ are gyromagnetic ratios of the nucleus 
and electron, S = ^ct and I is the nuclear spin (3/2 in 
this case). 

Our plan is to first express this operator in second- 
quantized form (in the electron basis) and then transform 
to the CF basis. 

Consider now the quantum well, which we take to be 
infinite in the x — y plane and of width w in the z- 
direction. In the absence any nuclear potential the single- 
particle wavefunctions will be given by 



*(k,s)(r,2,s) 



i^±{z)Xs 



(223) 
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where, Xs is a spinor with s = ±1/2, r, Hke k hes in the 
X — y plane, i!±{z) is a function hke •y/2/wsin(7rz/w). 
Note that in the z-direction there is essentiaUy just one 
wavefunction, other excitations being too high in energy 
to interest us. We now turn on the nuclear potential, 
whose effect will be to modify the states to 



*(k,s) (r, z, s) = £'''■'" ii(k^5) (r, z)xs 



(224) 



where, by Bloch's theorem, Uk,s(rj^) is invariant under 
translations by the unit cell in the x — y plane. The unit 
cell here extends into the z-direction from one end of the 
well to the other and may contain many nuclei. In the z- 
direction, Uk,s(r, z) will vary and possibly vanish rapidly 
beyond the ends of the well. We will assume U(k.s)(r, s) 
is normalized to unity over the sample. (This means it 
contains all normalization factors like 1/vL that we may 
usually associate with the plane wave.) The only thing 
of importance is that the states are still labeled by (k, s). 
As we shall see, detailed knowledge of U(k.s)(r, s) will not 
come into play. 

In second quantization we define a spinor field operator 



vl/(r,2) = ^d(k,s)e' 

s.k 



''"(k,s)(l',s)Xs 



(225) 



(Note that this operator is good for low energy physics 
in which higher excited states in the z-direction are ig- 
nored. We could add all these to the sum over states to 
obtain an operator that had Dirac delta function anti- 
commutation relations. However if wc limit ourselves to 
the low energy sector as defined above, this is a waste.) 

The spin density at the origin becomes 

(226) 



where 



"(k..)(O)-"(k,.)(0,0,0). 



(227) 



Now we make an assumption that is often made: 
W(k,s)(0) does not depend on s or k. We will simply 
call it u(0). Given this assumption 



s(o) = K0)P E 4',.'^rfks 



k,k'. 



k,qs,s' 

|«(0)p^5(q) 



(228) 

(229) 
(230) 



where 5(q) is familiar expression for spin density at mo- 
mentum q, the effect of the nuclear potential being en- 
coded in |w(0)p. Note that d]^ ^ still creates particles in 
states that are solutions to the nuclear and quantum well 



potential and not plane waves. However the commuta- 
tion rules of the d's are canonical and the spin density 
operator will obey the usual commutation rules 



[5"(q),5^q')] =i£'^''=5^(q + q') 



(231) 



When we go to the CF basis, we assume we will see 
the same change as in the case of the charge density: 



E^U'^k-E'^U'^ke^' 



(232) 



where I have included the factor e^'^ ' /^ so as to work 
with the projected density and not the magnetic trans- 
lation operator. Thus we write 



CF 



The 



S(0) 



factor 



^(0)P E 4^.V^'< 



k,k's,s' 
^ gJkxk'iV2g-|k-k'p/V4 



ikxk';2/2 -|k-k'p/^/4 



(233) 



which 



is just e^'^^^e~'^ ' /"* ensures that the above spin oper- 
ators (without the |u(0)p) have the right commutation 
relations ( Eqn. (28) of Moon et aE3) among themselves 
and the projected charge density. (The same criterion 
was used by Murthy in Ref.c3 in the small q limit.) 

We are now ready to eliminate |u(0)p in terms of a 
measurable quantity. The hyperfine interaction takes the 
form of an average field (B) acting on the nuclei: 






(234) 
(235) 



ks 



where I have used the fact that by symmetry, only k = k' 
and Sz can have mean values. 

Assume we are in a fully polarized state. Then the 
Knight shift is readily found to be 



K' 



;-fe-fnN\uiO)\' 



(236) 



where N is the total number of particles. For future use 
we invert this to write 



.(0)p = 



SisT'^ 



2-ieluN 



(237) 



b. Calculation of l/Ti 

Considp, now the relaxation rate. By the standard 
procedures one arrives at the following expression for 

l/Ti: 



1 _ lEmn^mn(£^m-£^n)' 

Ti 2 



Y.EI 



(238) 



W^n = 2^^^ |(TO^|7J,,/|n/)|2,5(£;f - Ei) (239) 



F I 
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where the bar indicates a thermal average, (m, n) label 
nuclear spin states, (/, F) denote many-body fermionic 
states, the energy difference. En — Em, between nuclear 
spin states has been neglected in the energy conserv- 
ing delta function, and H^f is the hyperfine interaction, 
which now know how to write in the CF basis. Thus we 
have 



7. Symbols 



w„ 



2tt 



87r7e7„|u(0)p 



X 



J2\J2{m\lM{F\SS^mi)\'SiEj. ^ Ej) (240) 



F,I a 

and we arrive at 



l/T, = ni 



/OO 



{0)S^'it)\I) (241) 



where S{t) is the Heisenberg operator at time t. 

Her e is a brief explanation of some steps leading to 
Eqn. ( ^4l[) . First I have used 



T E^ 



(242) 



where the transverse delta function S'^ means that a ^ z 
and then done some standard manipulations. 

To establish Eqn. (242), we need to invoke the follow- 
ing facts. 

• The kets \m) and \n) are eigenstates of the nuclear 
hamiltonian AI^ where A is some constant. 

• The factors {En — E,n)Ia or {En ~ E,n)Ia' can be 
traded for commutators of some other Ip or Ip' with 
A Iz- The result is zero if either a or a' equals z, 
which explains the transverse delta function that 
emerges. 

• The previous step allows one to invoke complete- 
ness and reduce the numerator to the trace of the 
hamiltonian {AIzY ^ which then cancels the denom- 
inator. 



Next we recall that 



k,s,k' 



s^'{t) = - ydi,(t)d,^M 



,ikxk';^/2 -|k-k'|^i^/2 



4.(i)=4.(0)e-*^=('^)* 

(and similarly for Sy^ ) in the HF approximation. In 
the above, £s(k) is the HF energy of a fermion of spin s 
and momentum k. 



To arrive at Eqns. ( |155H157| ) in the text, we do the 
integral over t (obtaining a delta function), use the ex- 
plicit expressions for spin operators, resort to standard 
HF factorization of the quartic operators and perform a 
standard change of variables in the measure, and finally 
eliminate |u(0)p via Eqn. (p37|). 



Symbol 


Significance 


ly 


= p/{2ps + l)=filling fraction 


P 


Number of CF LL's 


2s 


Number of vortices attached 


c2 


2ps/{2ps+l) 


B* 


Reduced field seen by CF = B/{2ps + 1) 


A* 


Reduced potential by CF = A/{2ps + 1) 


I 


electron magnetic length 


I* 


CF magnetic length = l/^/l — (? 


n* 


Velocity operator for CF 


p(q) 


Electron density in CF basis, 0{qV) 


x(q) 


Constraint in CF basis, 0{qi) 


p(q) 


Electron density in CF basis. 


x(q) 


Constraint in CF basis. 


P^(q) 


p(q) — c^/x(q) (preferred charge) 


/ 


Vortex form factor (e-9''V(8ps)) 


A 


Defined by v{q) = 2T:e^e~'i^^ /q 


b 


Defined by v{q) = 27re2e(«'^)' Erfc {qlh)/q 


^a,p 


Activation or polarizarion gap 


S 


A/(eVe/) 


1 / (2s) 


Defined by A^.p = eB*/{m^a,p 


^(2s) 


Defined by 1/mtp' = {eH/e)&^,;^ 


nor 
"^a,p 


'ma,p/{me\/B{T)) (normalized mass) 


rrif. 


Electron mass in free space 


P 


Polarization 


S 


Number of spin up minus down CF 


E{S) 


Ground state energy density 


9 


g -factor of CF, taken to be .44 


Ip-r,!") 


Many-body CF state with 




p — r spin LL's and r spin-down LL's. 


£±{k) 


Hartree Fock energy for up/down spin 


6 


Tilt angle 


Pnin2 


One-particle Matrix element of pf 




between LL ni and ^2 
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